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Abstract 



Incoherent photoproduction of pions on the deuteron in the A(1232) resonance re- 
gion is investigated in order to study the effect of nucleon-nucleon (NN) and pion- 
nucleon (nN) rescattering in the final state. The elementary 7^ — > irN production 
amplitude is taken in the effective Lagrangian approach and contains besides the 
standard pseudovector Born terms the resonance contribution from the A(1232) 
excitation. It yields for the elementary reaction a good agreement with the experi- 
mental data from MAMI, TRIUMF and TAPS. 

Pion photoproduction on the deuteron is dominated by the impulse approximation 
where the pion is produced on one of the nucleons neglecting all final state interac- 
tions. The comparison of the impulse approximation with the available experimental 
data shows a significant overestimation of the data. 

Therefore, the major point of concern of this thesis was the inclusion of rescattering 
effects in the final 7riViV system which we have limited to the leading order contri- 
butions of two-particle interactions in the NN- and 7riV-subsystems. As models for 
the relevant two-body interactions we have used separable approximations which fit 
the phase shift data for NN and -kN scattering. We found that the influence of 
NN- and 7riV-rescattering effects on total and differential cross sections is signifi- 
cant. Inclusion of such effects leads to a much improved agreement with the existing 
experimental data. 

Since several experiments to measure the spin asymmetry of the total photoab- 
sorption cross section, which determines the Gerasimov-Drell-Hearn (GDH) sum 
rule, are now being performed or planned at different laboratories around the world 
(MAMI, ELSA, LEGS, GRAAL and TJNAF), the contribution of incoherent pion 
photoproduction to the spin asymmetry for the deuteron is evaluated with inclusion 
of NN and itN rescattering contributions. The effect of final state rescattering on 
the spin asymmetry for the deuteron is found to be quite important and should be 
included in forthcoming theoretical studies. 
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ZUSAMMENFASSUNG 



In der vorliegenden Arbeit wurde die inkoharente Photoproduktion von 7r-Mesonen 
am Deuteron im Bereich der A(1232)-Resonanz unter Beriicksichtigung der Endzu- 
standswechselwirkung von Nukleon-Nukleon (NN) und Pion-Nukleon (7riV)-Ruck- 
streuung untersucht. Die verwendete elementare Amplitude der Photopionpro- 
duktion am freien Nukleon berucksichtigt Born-Terme in pseudovektorieller -rrN- 
Kopplung sowie den A(1232)-Resonanzbeitrag. Durch Anpassung der Modellpa- 
rameter erhalt man eine gute Ubereinstimmung mit den experimentellen Daten von 
MAMI, TRIUMF und TAPS. 

Beziiglich der inkoharenten Reaktion am Deuteron ergibt sich im Rahmen der Stofi- 
naherung qualitativ eine zufriedenstellende Ubereinstimmung mit dem Experiment 
fur geladene Pionen hinsichtlich der totalen und differentiellen Wirkungsquerschnitte. 
Es zeigt sich aber eine systematische Uberschatzung von etwa 30 Prozent fur neutrale 
Pionen. Daher ist es das Ziel dieser Arbeit, den EinfluB der bisher vernachlassigten 
A^A^-Endzustandswechselwirkung und der 7rA^-Riickstreuung auf die Pionphotopro- 
duktion zu untersuchen. 

Bei der Behandlung der Drei-Teilchen-Dynamik im Endzustand beschranken wir uns 
auf Beitrage von Zwei-Teilchen-Wechselwirkungen in den NN- bzw. 7rA^-Subsystem- 
en. Als Modell fur die NN- und 7rA^-Wechselwirkung haben wir separable Potentiale 
benutzt, wobei die freien Parameter durch einen Fit der Streuphasendaten bestimmt 
sind. Dabei zeigt es sich, daB der EinfluB der NN- und 7rA^-Riickstreueffekte auf die 
totalen und differentiellen Wirkungsquerschnitte signifikant ist. Die Einbeziehung 
solcher Effekte fuhrt zu einer verbesserten Beschreibung der experimentellen Daten. 

Als Beispiel einer Polarisationsobservablen haben wir die in die Gerasimov-Drell- 
Hearn-Summenregel (GDH) einflieBende Beitrage der inkoharente 7r-Produktion zur 
Spin-Asymmetrie des totalen Absorptionsquerschnitts am Deuteron unter Beriick- 
sichtigung der NN- und vrA^-Endzustandswechselwirkung untersucht. Wir haben ge- 
funden, dass der EinfluB der Endzustandswechselwirkungen fur die Spin-Asymmetrie 
und das GDH-Integral wichtig ist und in den zukiinftigen theoretischen Arbeiten 
berucksichtigt werden muB. 



Mainz, Mai 27, 2002 

Betreuer: Prof. Dr. H. Arenhovel 
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Chapter 1 

Introduction and Motivation 



1.1 Introduction 

The present work is concerned with photoproduction of 7r-meson on deuterium. The 
7T— meson or pion is the lightest of all the strongly interacting particles and plays a 
central role in the physics of strong interactions. It was proposed in 1935 by Yukawa 
PP as the carrier of the strong interaction between nucleons (proton, neutron) in 
analogy with the photon in the electromagnetic interaction. The existence of the 
7r-meson was experimentally confirmed in 1947 by Powell and his collaborators [2j 
and the pion is now known to exist in three charge states, 7r + , 7r _ and tt°, with 
masses of 139.6 MeV/c 2 for charged pions and 135.0 MeV/c 2 for the neutral pion. 

Photo- and electroproduction of pions on a proton have been studied thoroughly 
both theoretically and experimentally. Beyond the pion threshold energy almost 
all reactions on the nucleon, either initiated by photons, pions or other hadrons 
with energies of a few hundred MeV (in the center-of-mass) are dominated via the 
formation of the first excited state of the nucleon, the A(1232) resonance (J = 
3/2, 1 = 3/2), which almost exclusively decays into a pion plus a nucleon. 

During the last years, pseudoscalar meson production in electromagnetic reactions 
on light nuclei has become a very active field of research in medium-energy nuclear 
physics with respect to the study of hadron structure. For the following reasons the 
deuteron plays an outstanding role besides the free nucleon. The first one is that the 
deuteron is the simplest nucleus on whose structure we have abundant information 
and a reliable theoretical understanding, i.e. the structure of the deuteron is very well 
understood in comparison to heavier nuclei. Furthermore, the small binding energy 
of nucleons in the deuteron, which from the kinematical point of view provides the 
case of a nearly free neutron target, allows one to compare the contributions of 
its constituents to the electromagnetic and hadronic reactions to those from free 
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nucleons in order to estimate interaction effects. 

The basic interaction in photoproduction is as follows: a photon is incident on a 
target nucleus and interacts with its constituents. As a result, a pseudoscalar meson 
is produced along with other particles. Two kinds of processes depending on the 
nature of the other particles produced in this interaction are found: coherent and 
incoherent processes. 

In the coherent process, the meson is produced with the target nucleus maintaining 
its initial character. Thus, the interaction starts with a photon and some nucleus, 
and ends up with a meson and the same nucleus, i.e. 7X^4 — > tiXa, where A is the 
mass number of the target nucleus. The process is labeled "coherent" because all 
nucleons in the nucleus participate in the process coherently, leading to a coherent 
sum of the individual nucleon contributions. 

In the incoherent process, the nucleus ruptures and thus fails to maintain its initial 
identity. The meson is produced in association with a nucleon (or an excited state of 
the nucleon) and some new recoil "daughter" hadronic system. Thus, the interaction 
starts with a photon and some nucleus and ends up with a meson, a free nucleon (or 
an excited state of it) and a new hadronic system, i.e. ^Xa — > wNXa-i- The process 
is labeled as "incoherent" because it occurs in kinematic and physical circumstances 
similar to those of the process that produces a meson from a free nucleon. 

1.2 Review of Previous Work 

The electromagnetic production of pions on the free nucleon, including photopro- 
duction and electroproduction, has long been studied since the pioneering work of 
Chew, Goldberger, Low and Nambu (CGLN) As a result, an enormous amount 
of knowledge has been accumulated. Recently, theoretical interest in these reactions 
was revived by the new generation of high-intensity and high duty-cycle electron 
accelerators. With the developments of these new facilities, it is now possible to ob- 
tain accurate data for meson electromagnetic production, including spin-dependent 
observables. Extensive work during these more than fourty years (see for example 
Pj-|17j) indicates that, below 500 MeV incident photon energy, the mechanisms of 
the 'jN — > ttN reaction are dominated by the Born terms and the A(1232) excita- 
tion. 

In this work, we go a step further by studying photoproduction of pions on the 
deuteron. First investigations on photoproduction of pions on the deuteron go back 
to the early fifties [T%1 H§] with view on the general structure of spin flip and no 
spin flip amplitudes. Later, a more systematic calculation of pion photoproduction 
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on the deuteron was done by Laget [121 HB1 120] and Blomqvist and Laget [TT] . In 
their work the influence of pion rescattering and NN final state interaction is in- 
cluded within a diagramatic ansatz. They used the elementary photoproduction 
operator of Blomqvist and Laget [11 as input in their calculations. At the time of 
these calculations, a comparison with experimental data was possible only for 7r~ 
production since data for tt + and tt° production in the A(1232) resonance region 
were not available. The agreement of their predictions including final state inter- 
actions with the experimental data of the reaction jd — ► n~pp [21] is quite good. 
They found that the final state interaction effects are small for the charged pion 
photoproduction reactions in comparison to the neutral channel. In recent years, 
experimental data for tt° photoproduction on the deuteron have become available 
[22] • As mentioned in Ref . [22] > the predictions from [TT) IT2"1 IT3*1 |2*U] are significantly 
above the data. A possible reason for this may be that they used the Blomqvist and 
Laget parametrization JJJ of the elementary photoproduction amplitude which is 
not able to describe the neutral pion photoproduction from the proton. 

Pion photoproduction on the deuteron in the impulse approximation has been stud- 
ied by Schmidt et al. [23 3 neglecting all kinds of final state interactions and other 
two-body operators. They constructed an effective Lagrangian model for the process 
on the free nucleon and used it in their calculations on the deuteron. Since data for 
7r + and 7T° photoproduction on the deuteron were absent at the time of these calcu- 
lations, the authors could not compare their predictions with experimental data for 
these channels. A comparison with experimental data for ir~ production, showed 
a slight overestimation of the data. They reported that the reason for that is an 
overestimation of the elementary reaction on the neutron. A comparison between 
their predictions and the recently measured experimental data for 7r° photoproduc- 
tion on the deuteron [22] is now possible. One can see in [22] that this prediction 
can hardly provide a reasonable description of the data for this channel. As already 
noted in [T2"l [2H [2H], the effect of NN rescattering is important in the incoherent 
pion photoproduction on the deuteron, especially for 7r° production. 

Levchuk et al. [24J studied quasifree tt° photoproduction from the neutron via the 
d(7, n°)np reaction using the elementary photoproduction operator of Blomqvist 
and Laget JJJ. The contributions from the pole diagrams as well as one-loop dia- 
grams both with np and ttN rescattering were taken into account. They wanted to 
explore the possibility of measuring the i?i + /M 1+ ratio via photoproduction from 
quasifree production on the neutron. The isospin 1 = 3/2 component of this ratio 
characterizes the relative strength of the recently much discussed (see e.g. [2*H] l2*7j) 
quadrupole i?2-excitation of the A resonance. The idea was that the 71(7, ir )n reac- 
tion would be very useful for the isospin separation of the multipoles. In agreement 
with the results from the Laget model [12], Levchuk et al. find that the largest 
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effects disturbing the extraction of the multipoles for quasifree neutrons arise from 
the np final state interaction. They predict that these effects lead to a strong re- 
duction of the cross section at pion forward angles, but are much less important 
for backward angles. They found also that the correction due to np rescattering 
decreases with increasing pion angle and becomes to be less than 8% at 6 W > 90°. 
Furthermore, they pointed out that the contribution of the proton pole diagram 
and the one of ttN rescattering are negligible. The experimental data from Ref. [22] 
for the (i(7, 7T°)np reaction qualitatively support this prediction since the disagree- 
ment with the spectator approach is most severe at pion forward angles but less 
pronounced at backward angles. However, a comparison of the data to the Laget 
model including final state interactions shows some unexplained reduction of the 
cross section at backward angles. 

Recently and during the calculations of this work, Levchuk et al. [22] modified their 
theoretical predictions which have been done in (2^ using a more realistic version 
of the elementary production operator and including the charged pion production 
channels which are not included in their old calculations [21]. The elementary pro- 
duction operator is taken in on-shell form and calculated using the SAID [2E] and 
MAID [17 multipole analyses. The authors in [22] studied the semi-inclusive re- 
action g?(7, tt)NN in the A resonance region including pole diagrams and one- loop 
diagrams with NN and irN rescattering in the final state. Their predictions for total 
and differential cross sections including final state interactions show good agreement 
with the experimental data. However, Fig. 11.11 shows that a big difference between 
the theoretical predictions from [22] (dashed curve) and [21] (solid curve) is found 
for charged pion production reactions at pion forward angles in the impulse ap- 
proximation. Furthermore, the predictions from j2B] for the total cross sections for 
charged pion channels are found to be higher than the ones from [22] in the impulse 
approximation. 

1.3 Motivation 

The main motivation for studying photoproduction of pions on light nuclei is to 
obtain information on the elementary process on the neutron. The main goal of 
this thesis is to investigate incoherent pion photoproduction on the deuteron in the 
A(1232) resonance region in order to study the effect of NN and irN rescattering in 
the final state. We include besides the pure impulse approximation, the two-body 
t-matrices from NN and ttN rescattering in the final state. 

As already mentioned above, a big difference is found between the theoretical pre- 
dictions from [22] and [2B] in the impulse approximation. Therefore, it is interesting 
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Figure 1.1: Differential cross section for 7r~ photoproduction on the deuteron in the 
impulse approximation from [22] (solid curve) in comparison with the results from |25j 
(dashed curve). 



to check in this work where this big difference comes from. Since the authors in 
|2*B] pointed out that the main difference between their calculation and the one of 
Ref. is that a more realistic version of the elementary production operator is 
used, we will also examine in this work the use of different pion photoproduction 
operators. 

Furthermore, the comparison between the theoretical prediction from |2H] for the 
total cross section of 7T° photoproduction on the deuteron in the impulse approxi- 
mation and the TAPS data [22] (see Fig. ll.2|) gives a clear indication that the effects 
of final state interaction may be important. It was found that the final state in- 
teraction effects are significant in the case of coherent pion photoproduction on the 
deuteron (see for example [TU ESI EDI EI])- This means that one needs a reliable 
description for the rescattering process. 

Recently, several experiments to measure the spin asymmetry of the total photoab- 
sorption cross section, which determines the Gerasimov-Drell-Hearn (GDH) sum 
rule, on the proton and on the deuteron are now performed or planned at different 
laboratories around the world (MAMI, ELSA, LEGS, GRAAL and TJNAF). The 
A2 collaboration has prepared an experiment on the proton and on the deuteron 
at the Mainzer Microtron MAMI, that shall be completed for the higher photon 
energies at ELSA in Bonn. This makes the theoretical investigation of the spin 
asymmetry and the corresponding GDH integral particularly interesting. Therefore, 
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Figure 1.2: Total cross section for tt° photoproduction on the deuteron in the impulse 
approximation from [2B] as a function of the photon energy in the laboratory frame. Data 
points are from TAPS |22j . 



we investigate in this work the influence of NN and ttN rescattering effects on the 
spin asymmetry and GDH sum rule for the deuteron. 

Therefore, our aim is to construct a model for the reaction 'yd — > irNN in the 
A (1232) resonance region by incorporation the leading contributions from NN and 
irN final state interactions in time-ordered perturbation theory. 



1.4 Outline of this Thesis 

This thesis is organized as follows. In chapter we present the effective Lagrangian 
model of the elementary pion photoproduction process on the free nucleon which 
we use as input in our calculations on the deuteron. We explain how the reac- 
tion operator for this process is constructed. The nonresonant amplitudes and the 
contribution of the A (1232) resonance are given in an arbitrary frame of reference. 
We end this chapter with a discussion of our results for differential and total cross 
sections and compare with experimental data. 

Chapter |3] is devoted to the central topic of this work. The treatment of the 
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'yd — > ttNN amplitude, based on time-ordered perturbation theory with the in- 
clusion of iViV and ttN final state interactions is developed in this chapter. The 
transition matrix element is explicitly described with the inclusion of the NN- and 
7riV-rescattering. 

In chapter 0J we will present our main results together with a comparison with ex- 
perimental data and other theoretical predictions. We discuss in this chapter the 
effects of rescattering on total and differential cross sections. The contribution of 
incoherent single pion photoproduction to the spin asymmetry and the correspond- 
ing GDH sum rule for the deuteron is also discussed. Finally, we conclude and 
summarize our results in chapter Future considerations are also given in this 
chapter. 

For the convenience of the reader, eight appendices are given at the end of the thesis. 
Appendix El El and contain the general notations, formalism and useful formulas 
for the process on the nucleon. In appendix El we give the parametrization of the 
deuteron wave functions for the Bonn potential which we use in our calculations. 
In order to compare our results for differential cross section of 7r° photoproduction 
reaction on the deuteron with the experimental data from Ref. [221, transformation 
formulas from the laboratory frame of the deuteron to the 7./V center-of-mass frame 
are given in appendix El In appendix |Fl we study in detail the NN and ttN 
scattering matrices which we use as input in our predictions and also give a solution 
for the two-body scattering matrix using separable two-body interactions. The 
parameters of the form factors of these separable models are given in appendix IU1 
and El 
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Chapter 2 



Pion Photoproduction on the 
Nucleon 

Starting point of the construction of an operator for pion photoproduction in the 
two-nucleon space is the study of the elementary process, i.e. pion photoproduction 
on the free nucleon. This process is usually labeled as "elementary" to distinguish 
it from the same process from a nucleus. In the elementary process a photon is 
absorbed by a free nucleon (a proton or a neutron) to yield a 7r-meson in addition 
to a nucleon. 

The electromagnetic production of pions on the free nucleon, including photopro- 
duction and electroproduction, has long been studied since the pioneering work of 
Chew, Goldberger, Low and Nambu (CGLN) j3J. An enormous amount of knowledge 
has been accumulated as a result. Kroll and Rudermann [3] were the first to derive 
model-independent predictions in the threshold region, so-called low-energy theo- 
rems (LET), by applying gauge and Lorentz invariance to the reaction ^yN — > ttN. 
The general formalism for this process was developed by Chew et al. [3 (CGLN- 
amplitudes). Fubini et al. jS] extended the earlier predictions of LET by includ- 
ing also the hypothesis of a partially conserved axial current (PCAC). In this way 
they succeeded in describing the threshold amplitude as a power series in the ratio 
m = H 2 - up to terms of order m 2 . Berends et al. [6 J analysed the existing data in 
terms of a multipole decomposition and presented tables of the various multipole 
amplitudes constructing in the region up to excitation energies of 500 MeV. 

For more than twenty years, the standard model for pion photoproduction on the 
nucleon has been the model of Olsson and Osypowski E] , who emphasized the im- 
portance of Watson's theorem as a requirement to be obeyed by the electromagnetic 
multipoles below 7T7r-threshold. In practice, a model that has been more extensively 
used for comparison with data is the one of Blomqvist and Laget [TT] which is a 
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non-relativistic reduction of the model of Olsson and Osypowski. It has been con- 
structed in a general frame of reference, but Blomqvist and Laget used different A 
parametrizations for neutral and charged pion photoproduction. These parametriza- 
tions give a satisfactory fit to the amplitude for charged pion photoproduction, but 
it is not able to describe the neutral pion photoproduction from the proton. More 
than ten years ago, another model for pion photoproduction on the nucleon in the 
A region has been proposed by Nozawa, Blankleider and Lee This model has 
been extended to pion electroproduction and has also proved to be successful for 
photo- and electroproduction on the free nucleon. 

Garcilazo and Moya de Guerra ^E] have been constructed a model for photo- and 
electroproduction on the nucleon. This model is applicable from threshold through 
the first and second resonance regions. It contains the Born terms, the p and uj vector 
mesons, the A and the Roper resonances as well as the resonances Sn, D 13 , S31 and 
D 33 . They found a good agreement with experimental data up to 1 GeV. Recently, a 
unitary isobar model for pion photo- and electroproduction off the nucleon has been 
developed for nuclear applications at photon energies up to 1 GeV [Ej. This model 
contains Born terms, vector mesons and nucleon resonances and is constructed in the 
ttN center-of-mass frame. Within this model they have obtained good agreement 
with experimental data for pion photo- and electroproduction on the free nucleon. 

In this work we will examine the various observables for pion photoproduction on 
the free nucleon using the effective Lagrangian model of Schmidt et al. [23] which we 
will briefly outline in this chapter. The main advantage of this model is that it has 
been constructed in a general frame of reference and therefore can be applied directly 
to the electromagnetic photoproduction of 7r-mesons on nuclei. This model contains 
besides the standard pseudovector Born terms the resonance contribution from the 
A(1232) excitation. Kinematical and other useful formulas for pion photoproduction 
on the free nucleon, which diagramatically shown in Fig. 12.11 are given in appendix 

m 




Figure 2.1: The elementary process — ► ttN. A solid, dashed and wavy line represents 
a nucleon, pion and photon, respectively. 
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2.1 The Photoproduction Operator 

The T-operator describing transitions between asymptotically free states is given 
in terms of the interaction Hamiltonian Hint between all the involved particles as 
follows [S3 

T — 'Hint + Wint— t— — —'Hint , (2-1) 

E — TL + ie 
which can also be re-written as 

T = Hint + H int — — To — ; — ~1 > (2.2) 
-fc — Ho + «e 



where 7io is the free Hamiltonian. The on-shell matrix element T/j of Eq. ()2.2)1 is 
given in terms of the Hamilton operator Hint by 

Tfi = (f\T\i) = (f\H int \i) + J2(Wint\a) I . T ai , (2.3) 

a 

where the energy eigenvalues E a are given by 

H \a) = E a \a). (2.4) 

Making one-iteration approximation and keeping terms from the second order we 
obtain the following expression for the Tyj-matrix 

if = (f\H int \i)+Y,(f\ Hint \oc) E _l +ie (a\H int \i). (2.5) 

The complete Fock space of the system contains N-, tcN-, tctiN-, NNN-, ttNNN-, 
A-, Ti A-, etc states. 

For the process in our case, the form of the initial photon-nucleon state \i) and the 
final pion-nucleon state |/) is specified by the asymptotical states as follows 

\i) = \N,r,Pim t , ke) , (2.6) 
\f) = \N',7r;p 2 m t ,,qpL), (2.7) 

where p\, k and q are the momenta of initial and final nucleon, photon and 
meson, respectively. The isospin projection of the produced pion is given by /i, 
the polarization vector of the incoming photon by e* and rrit and are the isospin 
projection of the initial and final nucleon, respectively. The states of all particles are 
covariantly normalized (see appendix IA.2|) . The individual terms of the Tjj-matrix 
for pion photoproduction on the free nucleon are shown in Fig. 12.21 In the following 
we will evaluate the full interaction Hamiltonian of Eq. ()2.5|) in more details. 
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Figure 2.2: Diagrams for the elementary process 7iV — > ttN: (a) the Kroll-Rudermann 
graph, (b) and (c) the two time-ordered contributions to the direct and crossed nucleon 
pole graph, (d) and (e) the two time-ordered contributions to the pion pole graph, (f) and 
(g) the Z-graphs and (h) and (i) the A(1232) resonance graphs. Lines descriptions as in 
Fig.Eni 

2.2 Interaction Hamiltonian 

The general form of the interaction Hamiltonian of the involved particles, i.e. nu- 
cleon, A, pion and photon is described by the operator 7ii nt which is given by 

'Hint = 'Hem + ^vrV + T^TrjVA , (2.8) 

where TC em , 'H 7T n and TC^na are the Hamiltonians of the electromagnetic interaction, 
the 7riV-coupling and the coupling of the A resonance to the irN system, respec- 
tively. In the following we begin by evaluating the Hamiltonians which describe the 
interaction between pions, nucleons and photons as well as the contribution of the 
A(1232) resonance. 



2.2 Interaction Hamiltonian 
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2.2.1 The Electromagnetic Interaction 

The electromagnetic interaction 7i em contains in our case the coupling of the photon 
field to the free nucleon, pion and A fields 

Km = 7~L~fN + 7^77T + T^-ynN + 'H-yNA ■ (2-9) 

The Hamiltonian corresponding to the absorption of a photon at a nucleon or a pion 
(see Fig. 12. 3 J) are given, respectively, by 



Ti lN = -ej«(f)i(f)'7f(f) (2.10) 
and 

H 7 „ = ~ f d 3 x J2(-r^fiA(x) [(v$ M (£)) (x) - (v$_„ (£)) $ M (£) 

(2.11) 

where A(x), *S?(x) and $^(x) are the field operators of the photon, nucleon and 
pion, respectively. More information about these operators is given in appendix IU1 
7 represent the Dirac matrices and e denotes the elementary charge. 





Figure 2.3: The jN- and "fir- vertices. 

Using field quantization (see appendix EJ) the matrix elements corresponding to the 
two diagrams of Fig. 12.31 are given by 

(N';p 2 m' t \H yN \N,r,Pim t ,ke) = -(2n) 3 5 3 (p 2 - pi - k) 

x 2jk/~( m < I (e(p2+Pi)+i(e + K)ax kj -e \ m t ) (2.12) 

and 

{ir-q'yL^H^irAvM = -{2Tif5\q'-q- k)5^e/l(q + q') ■ e, (2.13) 

where is the nucleon mass and a are the Pauli spin matrices, e and k denote 
nucleon charge and anomalous part of the nucleon magnetic moment, respectively. 
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These are isospin operators of the nucleon and are given by 
e = |(l + r ) , 

k = ~ [ Kp (1 + T Q ) + K„ (1 - T )] , (2.14) 

where k p = |(k s + k v ) = 1.79 and n n = \{k s — k v ) = —1.91 are the anomalous 
magnetic moments of the proton and the neutron 1 in units of nuclear magnetons, 
respectively, and k s = —0.12 and k v = 3.70. The coupling constant ^ = 

In addition to the matrix elements of Eqs. (|2.12|) and (|2.13jl . the following matrix 
element of Ti. ni between a state of two-pion and one photon (see Fig. 12.4)1 is consid- 
ered 

(7r,7r';g/i,gV|^ 7 | 7 ;fce) = -(27r) 3 5 3 (<f + q - k)5^efx(q + q') ■ e . (2.15) 

This matrix element will be used later in the construction of the amplitude of di- 
agrams (d) and (e) in Fig. 12.21 The pionic current terms, which are given in Eqs. 
()2.13|1 and (|2.15)1 . contribute only to the photoproduction of charged pions. 



Figure 2.4: The 7r7-vertex. 



The 77riV Hamiltonian is given by 




(2.16) 



where is the pion mass and r represent the isospin matrices. We used the 

f 2 

irN coupling constant = 0.0735 which is given in Ref. [33 by fitting the rcN 
scattering data. The ^nN Hamiltonian is linear in photon and pion fields. This 
leads to a vertex, in which both photon and pion couple to the nucleon. The matrix 
element of T~l 77r N for the diagram in Fig. 12.51 is given by 2 

(N' \7r]p 2 m' t ,qfx\7C rKN \N,T,Pi'rrh,ke) = (2n) 3 5 3 (p 2 + q - pi - k) 

x — (mj | a - e [e,T+] | tth) . (2.17) 

^^The magnetic moments of the proton and neutron are /i p = 1 + n p and jtt n = K n , respectively. 
2 This term contributes only to the photoproduction of charged pions. 
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Figure 2.5: The jirN- vertex. 



Now we evaluate the fourth term in Eq. ([2.9)1 . In the description of the 7 iV A- vertex 
(see Fig one nas to take into account the magnetic dipole Ml and a possible 
electric quadrupole E2 excitation of the A resonance 




(2.18) 



Figure 2.6: The 7iVA-vertex. 



Since the strength of the electric quadrupole excitation E2 is much smaller than 
the magnetic dipole one (see for example (HOI EH] ) we will neglect it in this work. 
Following Weber and Arenhovel |35j . Wilhelm and Arenhovel [30J and Schmidt et 
al. > the jNA vertex reads 

ri-yNA = 7777 a AN ■ [kxe) T ANfi . (2.19) 

Here W^n denotes the invariant mass of the 7riV-subsystem and it is given by 
W nN = E N {q c . m ) +Un(q c .m.) , (2.20) 



where E N = a/ Mjy + q% m and io v = a/ m% + q 2 cm with the cm. pion momentum 
q c .m.- The transition spin (isospin) operator <jn A = <?atv (^na = t an ) is normalized 

as 

3 11 3 

(2 II van(t an ) II -) = -(-|| o-jva(tjva) II 2> = 2 - ( 2 - 21 ) 
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The energy dependent and complex coupling G^^W^n) is given as in jHHI by 
.mi , M/ \ _ / ^ M \W^ N )e^ M ^ w ^ for W nN > m7T + M N 



G aMWW-| dse , (2-22) 

where fi M1 (W n N) is given by 

^ M \W nN ) = /i + v* f — V + ^4 f — V (2.23) 

and the phase § M1 {W^ N ) by ggj 

$^(WW) = 9i 2 . (2.24) 

gA is the on-shell pion momentum in the irN c. m. frame on the top of the resonance, 
i.e., when the invariant mass W % n of the ttN state equals the mass of the A resonance 

W vN = luAqa) + E N (q A ) = M A . (2.25) 

It is given by 



\Wl N - ml - M%f - AmlM% 



The free parameters /io = 4.16, fi2 = 0.542, /i 4 = —0.0757, a± = 0.185 fm~ 3 and 
ci2 = 4.94 fm _1 are fitted to the experimental data for the M^ 2 -multipole of pion 
photoproduction [2E1 ED] • 

2.2.2 The ttN Interaction 

The irN interaction operator in pseudovector coupling is given by 

U*n = -— /"c/ 3 x^(f)7-7 5 r-^(f)V$(f). (2.27) 
m n J 

This operator is linear in the pion field operator Therefore, only one pion can 

be produced or absorped at the 7riV-vertex. Thus, only two possible diagrams as 
shown in Fig. 12.71 contribute. The evaluation of these two graphs yields the matrix 
elements 

(N' 1 ir;p2'm! t ,qfj,\7i vN \N;pimt) = -(27r) 3 5 3 (p 2 + q - Pi) 



x ^ (m' t | q- ar+ \ m t ) (2.28) 



2.2 Interaction Hamiltonian 
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Figure 2.7: The irN- vertices, 
for the emission of a pion and 

(N\7i]p2m' t ,qfi\H nN \N,7i 1 ,7i2-,pim t ,q 1 fi 1 ,q 2 fi2) = ~(^] 



6 2ifnN 



x [$ 3 (P2 - Qi -Pi)$ 3 (q- «2)<Wi a; fc(-) m ( m * I • | mt) 

+ 5 3 (p 2 - q2-pi)S 3 (q- gi)<^<fi(-rK I ^ • a | m t >] (2.29) 

for the absorption of a pion, where = \J m\ + qf is the energy of the pion with 
momentum Note, that in case of the pion absorption the vertex is evaluated for 
the transition of a two to a one pion state (see Fig. 12. 8 j) which will appear later in 
the construction of the matrix elements. 



N 

S 

S 



Figure 2.8: The 7riV-vertex in case of the pion absorption. 



2.2.3 The ttN A- Vertex 

Here we will evaluate the 7riVA-vertex which contribute to the amplitude of the A 
resonance. For the irN A- vertex which diagramatically is given in Fig. 12.91 we use 

mm nisi 

n nNA = -—F A (q 2 )(-Yr NA ^^a NA -q. (2.30) 



We have introduced a hadronic monopole form factor 

F A {q 2 ) = U NA ^-^. (2.31) 
Ai +q 2 
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Figure 2.9: The vrTVA-vertex. 

f 2 

The coupling constant = 1.393 and the cutoff A A = 315 MeV are fixed in Refs. 
|3*!H l3*7j to fit the irN scattering phase shift in the P33 channel and is also used in 
the calculations of this work. 

Now, using the Hamilton operator ()2.8)1 and taking into account all possible in- 
termediate states I a) in (|2.5jl . we can calculate the on-shell T^-matrix for pion 
photoproduction on the free nucleon by constructing the lowest order diagrams as 
shown in Fig. 



2.3 Construction of the Amplitude 

Using the previously mentioned electromagnetic and hadronic vertices, it is possible 
now to calculate the Tjj-matrix from Eq. ()2.5|) for photoproduction of pions on the 
free nucleon. Obviously, only the Kroll-Rudermann term 3 (T^ R ) contributes to the 
first term of the right-hand side of Eq. (|2.5jl . 

2.3.1 The Born Terms 

First, we consider the nonresonant amplitudes. These are referred to as the Born 
terms and they are dominant at low energy and for charged pion photoproduction 
still provide 50% of the cross section in the energy region of the A(1232) resonance. 

Using graphs (a) to (e) in Fig. 12.21 one finds in addition to the Kroll-Rudermann 
term (graph (a)) the direct and crossed nucleon pole terms (graphs (b) and (c)) and 
the two pion pole terms (graphs (d) and (e)). These terms are given, respectively, 
by 

Tjf = (2nf6"(p 2 + q-p 1 -k)^a-e[e,^], (2.32) 



7% = _(27r) 8 tf s (fc + 
J z 

3 Named also seagull or contact term. 
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r ^<? ■ q \2(p 2 + q) ■ e e + ia ■ k x e (e + k] 



( 2p 2 ■ e e + ia ■ k x e (e + k) ) r+c? • q 
+ ~ F fF F u) I ' (2 " 33) 



P2 + q - Pi - fc) 



T£ = (2 7 r) 3 5 3 (p 2 + g-p 1 -A;): 



x [ + ] [e, T+] , (2.34) 

where F^ = a/ Mjy + p 2 and cjp- = ^/ m 2 + p 2 are the energies of a nucleon and a 
pion with momentum p, respectively. 

To take the anti-nucleon terms in the propagators of direct and crossed nucleon pole 
graphs into account, we consider also the two Z-graphs given by diagrams (f) and 
(g) in Fig. 12.21 The matrix element using these Z-graphs (TS) is given by 

2% = (2n)H%p 2 + q-p 1 -k)^M N ^a-e 

E P2+q{ E P2+q + E P2 + Uq) E fo_%{E L + Fp> 2 - W 7 ) ^ 

We would like to note that for a better description of the real part of the M±+ 
multipole (see section 12.4)1 the suppression of the nonresonant background by the 
form factor is essential. Therefore, it is necessary to introduce a form factor 

Fs{q) = , (2.36) 

with the cutoff A# = 800 MeV to obtain a better description. 

2.3.2 The A (1232) Resonance Term 

The dominant non-Born contribution for photon energies up to 500 MeV is that 
of the P33 pion-nucleon resonance, the A(1232) resonance. Now, we will evaluate 
the contribution of the A(1232) resonance corresponding to the vertices given in 
Eqs. (12.19)1 and ()2.30)) . Using the nonrelativistic form of the A propagator, the 
various diagrams involving an intermediate A (1232) (see graphs (h) and (i) in Fig. 
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12. 2j) can be calculated. We obtain the following expression for the s and u channel 
contributions in center-of-mass frame 

Tn = (2n) 3 5 3 (p 2 + q-p 1 -k) 



x 



Fa (q 2 ) ef nNA G^(W wN ) {t% t } - |[rj, r ] g NA ■ qg AN -kxe 
m v 2JE~EJ 2 3 WW-M A + |r A (W, 



'ttN i 



| (0) ef vNA G%h(0) {t% t } + |[rt, r ] q AJV ■ fc x ea NA ■ q \ 
where Ma is the mass of the A resonance. 

The energy dependent width of the A resonance above pion threshold T A (W n N) is 
given by 

1 M N g 3 f 1 F 2 ( n 2 



T A (W N ) = < ^^W^W^aO?) for W 7rN >m 7T + M N 
* else 

In the u channel contribution which is given by the second part in Eq. (|2.37|) we 
take the values of the form factor Fa, the electromagnetic coupling an d the 
width of the resonance Ia at pion threshold. In the expressions of the hadronic form 
factors (Eq. (j2~3T|) ). the complex coupling Gfjj (Eq. RFM ) and the width of the 
A resonance (Eq. (|2.38|) ) we use the cm. pion momentum as given by the invariant 
mass of the 7riV-subsystem (see Eq. (12. 20(0 . 

In order to build the elementary pion photoproduction operator into the two-nucleon 
system, which will be necessary for the calculations on the two-nucleon space in the 
forthcoming chapters, we must construct the elementary amplitude in a general 
frame of reference. The Born terms given in section 12.3. II are already constructed in 
an arbitrary frame. In order to build the resonance term of Eq. 1)2.37)1 into the two- 
nucleon system we must re-write Eq. (j2.37j) in an arbitrary frame of reference. In 
Ref. f2Pj the following Galilei invariant transformations are found for the photon and 
pion momenta which have to be replaced by the relative photon-nucleon momentum 

k — >k lN (k,p l ) := — (2.39) 

Ma 

and respective pion-nucleon momentum 

— * { ~* ~* \ M N q-ujqp 2 fr)An . 
q — >qnN{q,P2) ■= -T7— • 2-40 

M N + ug 



Now, all contributions for the elementary photoproduction operator on the free 
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nucleoli, which in Fig. 12.21 are graphically represented, are calculated. The full Tfi- 
matrix is then given in a simple form by 4 

T wy = T™ + I% + Tf i + '[% + T£. (2.41) 



2.4 Results for Elementary Process 

In this section we will examine the various observables for pion photoproduction 
on the free nucleon using the effective Lagrangian model of Schmidt et al. [23 . 
This model contains besides the standard pseudovector Born terms the resonance 
contribution from the A(1232) excitation. The parameters of the A resonance are 
fixed by fitting the experimental data [21] for the M^ 2 multipole which is dominant 
in the A (1232) region. Due to the use of a constant A mass in the A propagator 
and a different cutoff Aa we had to increase from [HO] by a factor of 1.15 to fit 
the experimental M^ 2 multipole. 
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Figure 2.10: Real and imaginary parts of the M 1+ multipole in comparison with the 
calculation using the MAID program ^7]. The data points show the results of the VPI 
analysis [2B] (solution: September 2000) and the Mainz dispersion analysis (HDT) 27 . 



Fig. 12.101 shows our fit for the real and imaginary parts of the M l+ multipole in 
comparison with the MAID analysis ^7j, the Mainz dispersion analysis (HDT) [2*7j 

4 In this work no further contributions are used. This means in particular that terms from the 
lo- and p-meson exchange in the i-channel are not considered. 
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and the VPI analysis from the SAID program [2Bj- We see that the agreement of 
our results (solid curves) in comparison with data and theoretical calculation from 
MAID (dashed curves) is good. 



2.4.1 Differential Cross Section 

The cm. differential cross section for the transition from an initial photon-nucleon 
state \i) to a final pion-nucleon state |/) is given by 

_da_ _ 1 gM* y. a 

where u; 7 is the photon energy in the laboratory frame and m 7 = ±1. The magnetic 
quantum numbers of the target and the recoiling nucleons are respectively m s and 
m s i. The invariant ttN mass W n N is given in the cm. frame as 

W nN = E fl +u y , (2.43) 

and the absolute value of the pion momentum q is given by 



4 = vn»w.» y {2M) 



2W, 



ttN 



In Figs. I2.11H2.131 we compare our results (solid curves) for the differential cross 
sections with the MAID analysis (dashed curves) and with old and new experi- 
mental data from Mainz and Bonn. Most interesting for our analysis are the recent 
experimental data from the Mainz Microtron MAMI. Differential cross sections of 
both 7r + and 7r° photoproduction off the proton have been measured j2£l EH 139] . 
with high precision at all angles and photon energies for both channels. To get the 
full isospin decomposition we also had to include data for 7r~ photoproduction, for 
which we took the differential cross sections from j^U] ■ I n general, we obtain a good 
agreement with the experimental data for pion photoproduction on both the proton 
and neutron, especially in the energy range of the A(1232) region. 

For charged pion channels we see from Figs. 12.111 and 12. 121 a similar shape, but we 
find a somewhat different situations. Fig. 12.111 shows the result of our calculation of 
7r~ photoproduction on the neutron at six different values for the photon energy in 
comparison with the data from jlOj and also with the MAID analysis JE] . It can be 
seen that we fit this experiment quite well, in particular for a photon energy of 330 
MeV which is very near to the A(1232) region. Furthermore, the comparison with 
the MAID analysis at this energy is good. For higher and lower photon energies 
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Figure 2.11: Differential cross section for 772, — > pir~ at six different values of the photon 
energy in the laboratory frame. The solid curve shows the result of our calculations and 
the dotted one shows the results using the MAID program ^Jj. The experimental data 
are from Tokyo |40| . 



we note small discrepancies. In the case of high energies these discrepancies come 
from the fact that no other resonances besides the A(1232) are included in our 
calculation. It is also interesting to point out the importance of the Born terms in 
the charged pion photoproduction reactions in comparison to the contribution of the 
A resonance. These terms play an important role in the case of low photon energies. 

In Fig. I2.12l we compare our predictions for 7r + photoproduction on the proton with 
the more recent experimental data from the Mainz experiment j2Hj and the GDH 
experiment jSHj for different photon energies. We see that the agreement with both 
experiment and MAID analysis is good. We think that the agreement with 
experiment is in our case better than the MAID analysis for low and intermediate 
photon energies. For high energies small discrepancies are found since we assume a 
pure magnetic dipole transition in the 7^ A coupling. Moreover, the D i3 resonance 
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Figure 2.12: Differential cross section for "fp — ► nir + . Lines description as in Fig. 12.111 
The experimental data are from GDH [3H] and MAMI 



contributes with a non- vanishing term in this region. 

In case of the 7r° photoproduction on the proton, shown in Fig. 12.131 the situation is 
much more satisfactory. We see that the agreement of our calculation with the most 
recent experimental data from [SB] and [SHI is good. Small discrepancies between 
our calculation and the MAID analysis appear which very likely are due to the fact 
that no other resonances are included in our calculation. 



2.4.2 Total Cross Section 

The total cross sections for the different pion channels are shown in Fig. 12.141 and 
compared with experimental data. The total cross sections for 7p — > ir + n and 
jn — > n~p have a similar structure. In particular, the A peak is seen exactly at the 
same energy, 300 MeV. 



2.4 Results for Elementary Process 
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Figure 2.13: Differential cross section for jp — » pir°. Lines description as in Fig. 12.111 
The experimental data are from TAPS [HE] and GDH |39j . 



The comparison of our calculation with experimental data is carried out for photon 
energies up to 500 MeV. For higher photon energies the parameterization of the 

A resonance possesses no more valid. In general, we obtain a good agreement 

f 2 

with the recent experimental data using the small value -j^- = 0.069 for the pion- 
nucleon coupling constant. The agreement with the experimental data from [JUl 
and jH] for 7T~ photoproduction on the neutron is satisfactory. In case of the tt + 
photoproduction, the agreement with the recent experimental data from [39J is good 
up to photon energy of about 400 MeV. For high energies the Di% resonance, which 
is not included in our calculation, has non-vanishing contribution. The tt + data 
from jlHl are slightly overestimated in the resonance region by our calculation and 
also by the MAID analysis. 

At low photon energies, the charged pion reactions differ significantly in magnitude 
from the neutral pion reactions due to the Kroll-Rudermann term which does not 
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Figure 2.14: Total cross sections for pion photoproduction on the nucleon as a function 
of photon energy for the four physical reactions. The solid curve shows the result of 
our calculation, the dashed curve shows the results using the MAID program |17j and 
the dotted one shows the results of the VPI analysis of the SAID program [2B]- The 
experimental data are from GDH [U], Tokyo 001, TRIUMF gT] and TAPS 02]. 



contribute to the neutral pion reactions. The agreement of our calculation with the 
experimental data from [12] and [SHI hi case of the n° photoproduction on the proton 
is good and give a clear indication that the model which we used in our predictions 
can be applied directly to calculate the electromagnetic photoproduction of pions 
on the deuteron which we will study in the forthcoming chapters. 



Chapter 3 

Photoproduction of 7r-Mesons on 
the Deuteron 



In the previous chapter, we have constructed the elementary photoproduction op- 
erator on the free nucleon in a general frame of reference to use it as input in our 
calculations during the present chapter. This chapter is concerned with incoherent 
pion photoproduction on the deuteron which can be produced through the following 
three reactions 

jd — > ppix~ , 7<i — ► nmT + , 7c? — ► pun . (3-1) 



Let (u)y,k), (E?,d), (u}^,q), (Ep iy pi) and (Ep 2 ,p 2 ) be the four-momenta of the incoming 
photon, deuteron, the outgoing pion and two nucleons, respectively. The general 
expression for the unpolarized differential cross section of pion photoproduction 
reaction on the deuteron is given using the conventions of Bjorken and Drell [43] by 

1 fc\ \ — 5^4 / 1, 1 j n 1 ld 3 qd 3 Pl d 3 p 2 M 2 N 

da = (27r) [k + d — pi — p 2 — q) 



|uy — Vd\ 2 2u)tf Efi ± Ep 2 2ujy2Ej 

X 5 \^fm a ,m-ym d \ 2 y (3-2) 

smst,m~/md 

where m 7 is the photon polarization, the spin projection of the deuteron, s the 
total spin of the two outgoing nucleons, m s its spin projection, t the total isospin 
of the two outgoing nucleons, m t (= — /x) its isospin projection and w 7 and Vd are 
the velocities of the photon and the deuteron, respectively. \i denotes the pion 
charge. The states of all particles are covariantly normalized (see appendix IA.2|) . 
The reaction amplitude is given by 

•^sm s ,m^m d = (fiA, sm s , t, -fi\ M(q, fi, k, m 7 ) \d, m d ) . (3.3) 
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Choosing the z-axis in the direction of the incoming photon and isolating the az- 
imuthal dependence of the direction of pion momentum, we obtain the following 
general form for the reaction matrix 



Using parity conservation one can show that the .M-matrix elements obey the fol- 
lowing symmetry relation 



This symmetry relation reduces the number of complex amplitudes from 24 to 12 
independent ones. For their determination one needs 23 real observables since a 
overall phase remains arbitrary. 

The unpolarized differential cross section in Eq. (j3.2j) is a product of two terms. The 
first term is a kinematical one describing the invariant phase space. The second one 
describes the dynamics of the process and is determined by the transition matrix 
element M.. In the following we disscuss these two terms in more detail. 



Calculations of total and differential cross sections are carried out in the rest frame 
of the deuteron (laboratory frame) which is shown diagramatically in Fig. 13.11 In 
this frame the relative velocity of both particles in the initial state is unity in units 
of c. The coordinate system is chosen to have a right-hand orientation with z-axis 
along the momentum k of the incoming photon and y-axis parallel to k x q. The 
x-z-plane is then given by the outgoing pion, i.e. (fi n = 0, where 4>n is azimuthal 
angle of the pion. 

The four-momenta of the incoming photon and deuteron are given in the laboratory 
frame by fc M =(u; 7 , k) and <i M =(Md,0), respectively, where M d is the deuteron mass 




and uj 1 is the photon energy in the laboratory frame. Energy-momentum conserva- 
tion gives 




(3.4) 




(3.5) 



3.1 Kinematics 




k 



q + Pi+p 2 



(3.6) 



and 



uj^ + M d 



— ujf+ + E^ 2 — E ld , 



(3.7) 



3.1 Kinematics 
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Figure 3.1: Kinematics in the laboratory frame for the reaction jd — > ttNN. 

where E^ d is the total energy of the system. 
The invariant mass of the 7c?-system is given by 



The final state is determined by the four momenta of the two outgoing nucleons and 
the outgoing pion. Energies of these particles are then fixed by 



Finally, four-momentum conservation reduces the number of the independent vari- 
ables to five, out of nine. In this work, the pion momentum q, its angles 8 n and <f) n , 
the polar angle 8 Pr and the azimuthal angle <p Pr of the relative momentum p r of the 
two outgoing nucleons are selected as independent variables. We prefer this choice 
of variables, because in this case the kinematical factor, i.e. the phase space factor, 
does not have any singularities on the boundary of the available phase space, when 
p r —>■ (see also This means in particular that with another selection of the 







\/m\ + q 



(3.9) 
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free kinematical variables singularities in the phase space, i.e. kinematical factor, at 
the boundary of the physical region can occur (see for example |44|). 

The total and relative momenta of the final TVTV-system are given respectively by 

P = Pi+P2 (3.10) 

and 

Pr = g (Pi -n) ■ (3.H) 
The inverted relations read 

Pi = \ P + Pr (3-12) 

p 2 = \P-Pr. (3.13) 
The absolute value of the relative momentum p r is given by 



l E 2 NN (W 2 NN -AM 2 N 
E% N -P* cos* e PpT 



n — I NN V NN ^ 1V1 N) /o -,a\ 



where 



E NN = u-y + M d — = + Ep 2 , 

W 2 NN = E 2 NN -P 2 (3.15) 

and 9p Pr is the angle between P and p r . The absolute value of the total momentum 
P=\P\ is given from Eqs. (|3.fij) and ()3.10j) by 

P = | k-q \ . (3.16) 

The main features of the processes (|3.1j) will be investigated by considering the 
partially integrated differential cross sections d 3 a/(dfl n dq) and d 2 a/dQ n , which are 
obtained from the fully exclusive cross section 

d5(T =P± V \M^ I 2 C3 17) 

dQ p dtt n dq 6 ^ ' sm s ,m^m d \ \ ■ ) 

by appropriate integration. The phase space factor p s depends on the selection of 
the five independent variables in Eq. ()3.2|) . It is expressed in terms of relative and 
total momenta of the two final nucleons as follows 

1 r> 2 M 2 a 2 

(2tt) 5 \E a ( Pr + \P cos 6 Ppr ) + E fl ( Pr - \P COS 9p Pr )\ l^MdUq ' 1 ' ' 



3.2 Two-Nucleon Wave Function 
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whereas q between and q max - Since the limit value of both the numerator and the 
denominator are zero when p T — > 0, we can apply the L'Hospital's rule to obtain the 
exact limit value for the phase space as follows 

hm p s = Iim 



S Pr-o-£-\E f2 (p r + ±Pcos9p Pr )+E fl (p r -±Pcos9 Ppr )\ 

1 MW 
X (2vr) 5 16w 7 M d a;g- 
= 0. (3.19) 

The existence of this limit value ensures that the phase space factor has no singu- 
larities on the boundary of the physical region. 



3.2 Two-Nucleon Wave Function 

Using a covariant normalization (see appendix IA.2|) the deuteron wave function in 
the momentum space is given in the rest frame of the deuteron by 

E 

N m B 

where 



^m d mM = (2vr)^ V2Md Yl i L UL(p)Y LmL (p)(Lm L lm s \lm d ) . (3.21) 



L=0,2 



For the radial wave functions ui(p) we use the parametrization of the Bonn potential 
(full model) [15] (see appendix IDf 

-E^tH' (3-22) 
2 ^ Di 



V 7r *—f p 1 + mf 

i=i ' 

In case of the impulse approximation (see section l3.3.1j) . we use for the NN final 
state a complete antisymmetric NN plane wave. For the total spin and isospin 
part of the two nucleon wave functions we use a coupled basis (\sm s )\tm t )). The 
corresponding antisymmetric final NN wave function can formally be written as 

\p x p 2 ,sm s ,tm t ) = -L {\pi) {l) \p 2 ) {2) ~ {-) s+t \P2) {1) \pi) {2) ) \sm s )\tm t ) , (3.24) 
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where \f>i) describes a state in which a nucleoli "j" has the momentum pj. Only 
the t = 1 channel contributes in the case of charged pions whereas for ir° production 
both t = and t = 1 channels have to be taken into account. 



3.3 Matrix Elements for -> ttNN 

The transition matrix elements M.sm s ,m^m d are calculated in the frame of time- 
ordered perturbation theory, using the elementary operator introduced in chapter 
and including NN and irN rescattering in the final state. Scattering reactions are 
described in terms of probability amplitudes relating the initial and final asymptotic 
states of the combined system of projectile and target. In this section we consider 
the transition matrix element M. of Eq. ()3.2|) which describes the dynamics of the 
processes (13. lj) . These processes are diagramatically shown in Fig. 13.21 There, all 
contributions of two-body currents are neglected. All possible rescattering diagrams 
are shown separately in Fig. 13.31 The .M-matrix has the general form 

M = (irNN; qwp 2 ,a\e^(Q)\7d-,kd,0) , (3.25) 

where J M (0) is the current operator, a and (3 stand for the quantum numbers of the 
states asymptotically. 



N 



Figure 3.2: Diagramatic representation of the 7c? — > irNN amplitude including rescatter- 
ing in the final state and neglecting all contributions of two-body currents. 





T 




T 

71NN 








d 













We include in this work besides the pure impulse approximation (IA), the driving 
terms from NN- and 7rA-rescattering, so that the full transition matrix element 
reads 

A/fM = A/f(*M) IA _|_ A^M NN 1 Kj((tii) ttN (o 0(\) 

in an obvious notation. 

Now, we will consider successively all the three terms of Eq. (|3.26J) . i.e., the impulse 



3.3 Matrix Elements for — > ttNN 
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Figure 3.3: Diagramatic representation of final state interactions of the reaction 7c/ — > 
ttNN: (a) and (b) are "driving terms" from NN and irN rescattering in the final state, 
respectively. 



approximation, the NN final state interaction and the irN rescattering in more 
detail. 

3.3.1 The Impulse Approximation (IA) 

In order to qualitatively explain the approximations concerning the rescattering 
terms which are of important interest in this work and discussed below, we would 
like to demonstrate here some features of the reaction amplitude keeping in ()3.26j) 
only the IA term. In case of the 7c? — > ttNN reaction, the impulse approximation 
leads to the so-called spectator nucleon model. The Feynman diagrams of this model 
are shown in Fig. 13.41 The corresponding transition matrix M. IA is then given by 

M ia = T$1® + iVtV), (3.27) 

where the upper index indicates on which of the nucleons the elementary production 
operator acts. The operator M. IA contains contributions of pure single nucleon 
terms. 

Using the initial and final states from section 13.21 and the operator of the pion 
photoproduction in the two-nucleon space from Eq. ()3.27|) we can write now the 
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Figure 3.4: Feynman diagrams for the reaction 'yd — > ttNN in the impulse approximation. 

.M /j4 -matrix element using a covariant normalization for all particles (see appendix 
Z2I) as follows 



(pip 2 , sm s ,t, -fi\M IA (k, m 7 , q, p)\d, m d ) 



1 f d 3 p[ f d 3 p' 2 M N M A 



N 



X 



2 J (2tt)3 J (2tt)3 E n E p , 
^ (piP2,sm s ,t, -/i|.M /j4 (fc,m 7 , q,fj)\p{p2, lm s /,00) 



x (p^,lm s /, 00|d,m d ). (3.28) 
In the laboratory frame one finds for the matrix element ()3.28|) the following expres- 



sion 



(piP2,sm s ,t,-p,\M (k,m y ,q, p)\d = 0,m d , 00) 
1 



^< sms ,t,-//| [«ft|T«| -p 2 ) - (_)-+^|T®| -£» ^ ms/md (p 2 ) 

+ ((ftlTWl-^-H^lTWl-ft))^^)] |lm,,00>. (3.29) 

Note that the upper index is maintained on T nl . The reason for that lies in the 
fact that the T^-matrix contains spin and isospin operators, which still act on the 
spinors of the NN system. Using the symmetry properties 

{sm s \\ {1) a^\s'm sl ) = (-) s - s \sm s \a^ l (2) |s'm s ,) , (3.30) 

which applies also likewise to the r-matrices, the matrix element of Eq. 1)3.29)1 is 
given by the following expression 

Mty$_™„, = (p 1 p 2 ,sm s ,t,-{i\M IA (k,m y ,q,/j,)\d = 0,m d ,00) 



s/2j2(srn s ,t,-v\((pi\T$\ - p 2 )*m s „m d (P2) 
- (-) s+ *<^2lr«| -p 1 ># ms , ;md (pi))|l^,00) . 



(3.31; 
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The evaluation of the spin and isospin operators becomes very simple, because they 
act now on the nucleon "1" only. Using the Wigner-Eckart-theorem (see for example 
[IE]) one finds 

(sm s \al l) \lm s/ ) = (-) 5 " 1 3V2 (lvlm s ,\sm s ) j \ \ | } , (3.32) 

(tm t \rl^\00) = (-T5 til 8 mt ^. (3.33) 
3.3.2 The NN Final State Interaction 

Now we will consider in addition iViV-rescattering in the final state. The corre- 
sponding matrix element is given by the second term in Eq. ()3.2(i|) . The Feynman 
diagram for this process is shown in Fig. 13.51 for the case when the production op- 
erator acts on nucleon "1". Note that the second diagram when the operator acts 
on nucleon "2" is also taken into account in the calculation. 

The transition matrix element with iViV-rescattering of the diagram in Fig. 13.51 is 
given by 




P' 

Figure 3.5: Feynman diagram for the reaction jd —* irNN including iViV-rescattering in 
the final state. 

M^l^ md = (q,p^p 2 ;^,sm s ,t,-fi\M NN (k,m 1 )\d=0,m d ,00) 
1 x ^ [[[ d 3 q' d 3 p{ d 3 pj 
2 ^\^J" 2u q ,E pl /M N E p ,/M N 
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x (q,PiP2] V, sm a ,t, -fi\K NN (W NN )\q',p[p2] fi',s'm s ,,t', -//) 
x g^{E ld ,q,p'i,pL) 

x (q',p[p2^\s'm sl ,t',-fi'\M IA {k,m 1 )\d = 0,m d , 00). (3.34) 

Here TZnn denotes the half-off-shell iVTV-scattering matrix, Gonn (E^d, q, p[, P2) * s 
the free irNN propagator and .M /A -matrix is given in Eq. ()3.31|) . 

3.3.2.1 The Half-Off-Shell iViV-Scattering Matrix 

The NN dynamics in the final state is determined by the half-off-shell iViV-scattering 
amplitude TZnn- In the presence of a spectator meson it is given by introducing the 
relative and total momenta p r and P, respectively, of the two outgoing nucleons (see 
Eqs. (I3.10|) and (|3.11Jl ) and using the basis states 

{p r r (f) r \p r ,sm s ,tmt) = (sm s £m e \ JMj) ufj(j) r ) Y imt (j) r ) \£sJMj) \trrit) 

JMjlmi 

(3.35) 

by 



{q : PiP2]^,sm s ,t,-^\n N N{W N N)\q\p[P2]^\s'm s ^ = 2(2ir) <V/A'A't 

x 5\q' - q)5\P ! - P) J^^^ D s Z ms ,(W N NiPrX) , (3-36) 

with 

V S ^ s m s X W NN,Pr,P'r) = (Pr, Sm s , t - ll\K NN {W NN )\p' r , S77ly, t - (l) 

= J2^i£'m a m 3 ,(Pr,Pr) {W N N , Pr , P'r) , (3.37) 
JW 

where I is the orbital angular momentum and J is the total angular momentum of 
the two-nucleon system. The initial and final relative momenta of the two nucleons 
in the iViV-subsystem are given, respectively, by 
1 

Pr = -^{Pl-PL), 

Pr = ~ (Pi - #2)- ( 3 - 38 ) 

where p{ — p' r + — q) and p' 2 = —p' r + \ {k — q) with the integration variable p' r . 
The function Tf§ m#f?v (p n p' r ) is given by 

^ee'em s m s ,{Pr,Pr) = Yl sm s \J M j)(£' m e sm sl \J M j)Yl mt (p r )Y llmi ,(p' r ) ■ 

(3.39) 
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The half-off-shell partial wave amplitudes T u f ^{WnniPtiPI) were found for partial 
waves with total angular momentum J < 3, which are important in the energy 
region under consideration, by numerical solution of the Lippmann-Schwinger (LS) 
equation |Hj . In appendix |F] we give a detailed solution of the LS equation for a 
given NN potential model. In the calculations presented here, these amplitudes 
are obtained from the separable representation of the Paris NN potential |4*£| l4*H] , 
which reproduces its on-shell as well as off-shell properties. 



3.3.2.2 The nNN Propagator 

The free nNN propagator is given by 

= Ejd _^ {ql _ EN ; {pl) _ ENM) + M - (3-40) 
In the nonrelativistic limit it is given by 

g 0NN {E ld ,q, Pl ,p 2 ) = — -, (3.41) 



where p' is given by 



p' 2 = M N [E ld -u^-2M N - i ^^ \ . (3.42) 



The magnitude of the relative on-shell momentum of the two nucleons is given by 



Po = VM N (W NN -2M N ). (3.43) 



Collecting the various pieces and substituting Eqs. (|3.3fij) and (|3.41j) in Eq. (|3.34j) . 
we obtain the following expression for the transition matrix element with NN rescat- 
tering 

Mt\£%{KlPUP2) = E / d "Pr\^¥ V% ms ,(W NN ,Pr,ti) 

m s , J V Pi P2 

x M^\^ md {k,q^p' 2 ) , (3.44) 

p' —pl+ie 

where T>^ m ,(WNN,Pr,Pr) is given by Eq. f)3.37j) . The three dimensional integral 
over p' r is carried out numerically. 
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3.3.3 The ttN Final State Interaction 

In addition to the NN final state interaction we consider also in this work the ttN- 
rescattering in the final state. The corresponding matrix element in this case is 
given by the third term in Eq. ()3.26j) . The Feynman diagram for this process is 
shown in Fig. 13.61 for the case when the production operator acts only on nucleon 
"1". The second diagram when the operator acts on nucleon "2" is also taken into 
account in the calculation. 

The transition matrix element with 7riV-rescattering in the final state (see Fig. 13. 6|) 
is given by 




P2 

Figure 3.6: Feynman diagram for the reaction — * ttNN including 7T-/V-rescattering in 
the final state. 

M^iZ% d = (q,pm;fM,sm s ,t,-v\M* N (k,m J )\d=0,m d ,00) 
1 f f f d 3 q' d 3 p{ d z p' 2 

t^WfJJJ ^>E pi /M N E p ,/M N 
x {(q,PiP2;H,sm s ,t, -fx\n wN (W wN (p2))\q\p[P2;^',s'm s >,t', -//) 

G£S? 

x (q',p[p2^\s'm s/ ,t\-fi'\M IA {k,m 1 )\d = 0,m d ,00), (3.45) 
where IZ^n denotes the half-off-shell 7riV-scattering matrix at the invariant mass 

W« N {p 2 ) = ^(Efr+ujtf-fa + q)* (3.46) 

of the 7iN subsystem, (E ld) (t,pi,p£) is the free nNN propagator (see Eq. 

(jODjl ) and A^ /A -matrix is given in Eq. (}33Tj) . 
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3.3.3.1 The Half-Off-Shell vriV-Scattering Matrix 

Using the coupled basis states 

\q,p~[p 2 ;fi,sm s ,t,-fi) = ^ (^rn Sl ^m S2 \ sm s )(-m tl ^m t2 \t - fi) 

m si m tl ~ " " 

rn 32 rnt 2 

x \q,pi,P2-,lfi, |m Sl ,|m tl , \m S2 , \m t2 ) , (3.47) 

the half-off-shell 7riV-scattering amplitude H^n is given in the presence of a spectator 
nucleon by 

E- 

(q,piP2;t*,sm s ,t, -fi\n vN (W vN (p 2 ))\q\p[p! 2 ]fi\ s'm s ,,t', -fx') = (2tt) 3 -^ 

M N 

x m - Pi) E A Z>>' ( m s2,m t2 , m^miO 

m so mt mi m.i 

x (q,p 2 ; lfi, \m S2 , \m t2 \TZ vN {W vN (p 2 )) \q' , p 2 , Ifi', \m s ^ \m t ^) , 

(3.48) 

where m S2 {m s ^) and m t2 (m t ^) are the spin and isospin projections of the final 
(initial) nucleon which is interacting with the pion in the irN subsystem, respectively. 
The function A is given by 



,sm s s'm s , i \ \ ^ r r 



.1 1 ..11. 

x (\ m s[ \ m s' 2 I ^m t ' 2 I t' - fi) , (3.49) 

where m si and m tl are the spin and isospin projections of the spectator nucleon, 
respectively. 

Coupling the pion and nucleon isospins into the total isospin t of the pion-nucleon 
pair according to 

%p 2 ;\fi,\m S2 i\m t2 ) = ^{lfi]^n t2 | tm t ) \q,p 2 ;trh t , \m S2 , \m t2 ) (3.50) 

trh t 

and introducing the relative and total momenta of the final (initial) pion-nucleon 
subsystem which are given, respectively, by 

~ M N q- (=>, m nQ' - m w p^ \ 

Pr ~ M N + m w \ Pr ~ M N + m7V J 
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and 

P = q + p 2 (p> = q' + p^j (3.52) 
we can re-write Eq. ()3.48j) in the form 

E- 

{q,fiP2,^ sm s ,t, -ii\n nN (W nN (p 2 ))\q\p[p^fi',s'm s ^t',-fi') = 2(2%) 9 "'" 



M, 



N 



1V1 N 

E,sm 3 s'm , , \ 
{m S2 ,m t2 ,m s > 2 ,m t > 2 ) 



X 

J—J w 

irht t'm t i 

x (p r , ^m S2 Jm t \TZ 7rN (W 7rN (p 2 ))\p^ \m s ^t'm t ,) . (3.53) 
The function B is given by 

mm 1 1 

B TJii' ( m *2' m *2) = 2 m ' 2 I ^ ^ 2 m *' 2 ' ' ( ^ 3 ' 54 ' ) 

Using the partial wave expansion 

(p r 6<p\p r ,±m S2 ,tm t ) = ^ {]^rn S2 im l \ JMj) u\j{p r ) 

JMjimi 

x Yyjir) \i\JMj) \tm t ), (3.55) 

where £, J and Mj denote the relative orbital angular momentum, the total angular 
momentum and its z-axis projection of the ttN system, respectively, we obtain 

(Pr, ^m S2} tm t \n 7TN (W 7rN (p 2 ))\p^ \m s ^t'm v ) = QZ t m t X w ^N{p2),Pr,p') 
= S t'^m timt Sir i S me , me 5j,j5 £ljiA - Ij (^m S2 im e \ 

JMjIfht 

(3.56) 

Inserting Eq. ()3.56|) in (j3.53|) . we obtain the final form for the half-off-shell ttN- 
scattering amplitude. In the calculations presented here, the half-off-shell partial 
wave amplitudes ^^(W^-at^),^,^) were obtained for all 5*-, P- and D- waves 
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channels by solving the Lippmann-Schwinger equation [Uj for a separable energy- 
dependent nN potential built in Ref. [T^. This potential describes precisely the 
on-shell as well as the off-shell properties of the ttN scattering amplitude. More 
details about this potential is given in appendix |FJ 

Collecting the various pieces, we obtain the following expression for the amplitude 
1)3.45)1 with 7riV-rescattering in the final state 




x 




(3.57) 



where p' r {p) and p r are given in Eq. 1)3.51)1 and 



q'(p) = P-Pi 
#|(p) = -P- 



with pi 



Pi - k 



(3.58) 



Similarly to the case of iViV scattering, the three dimensional integral over p is 
carried out numerically. 
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Chapter 4 

Results and Discussion 



In this chapter we will present and discuss our results for pion photoproduction on 
the deuteron in the A(1232) resonance region including iViV and tcN rescattering in 
the final state. We will discuss the important observables that have been measured 
experimentally, such as the total and differential cross sections. Since several exper- 
iments to test the Gerasimov-Drell-Hearn (GDH) sum rule both on the proton and 
on the neutron are now in preparation or planned at different laboratories around 
the world (MAMI, ELSA, LEGS, GRAAL and TJNAF), the theoretical investiga- 
tion of the spin asymmetry and the corresponding GDH integral are nowadays a 
very interesting subject of research. This sum rule is derived under very general as- 
sumptions (Lorentz and gauge invariance, causality, relativity and unitarity). This 
makes its verification to be an important check of our understanding of the hadronic 
spin structure. The inclusion of final state interaction effects may play an important 
role on the spin asymmetry and the corresponding GDH sum rule for the deuteron. 
Therefore, we will also discuss in this chapter the influence of rescattering effects on 
the spin asymmetry and the GDH sum rule for the deuteron and, for comparison, 
for the nucleon. 

The discussion in this chapter is divided into three parts. In the first part, we 
present our results for the total cross section of the processes (13.1)1 in comparison 
with experimental data and other theoretical predictions. In the second part, we 
discuss the semi-exclusive differential cross section d 2 a/dQ 7T . We also compare our 
results for the differential cross section with experimental data and other theories. 
In the third part, we study the influence of NN and ttN rescattering on the spin 
asymmetry and the corresponding GDH sum rule for the deuteron. We would like 
to mention that the results presented here are calculated using the deuteron wave 
function of the Bonn potential (full model) [13] (see section l3~2*j) . 

To calculate the semi-exclusive differential cross section d 2 o / 'dQ n , the fully exclusive 
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differential cross section d b a / {dVL^dq^dVL^) (see Eq. 13.17)1 ) was integrated over the 
pion momentum q n and the polar angle 9 Pr and the azimuthal angle cj) Pr of the relative 
momentum p r of the two outgoing nucleons. The total cross section was calculated 
by integrating Eq. ()3.17jl over all the remaining independent variables, i.e., the pion 
momentum q n , its angles 6^ and 0^, the polar angle 6 Pr and the azimuthal angle 
4> Pr of the relative momentum of the two outgoing nucleons. These integrations are 
carried out numerically. The number of integration points was being increased until 
the accuracy of calculated observables becomes good to 1%. 



4.1 Total Cross Section 

We start the discussion with presenting our results for the total cross sections in 
Fig. 14.11 including NN and ttN final state interaction effects. 



300 



300 




_ D 



150 250 350 450 
c^ ab [MeV] 



150 250 350 450 
c^ ab [MeV] 




150 250 350 450 
c^ ab [MeV] 



Figure 4.1: Total cross sections for *yd — > ttNN reaction obtained within the impulse 
approximation (dotted curves) in comparison with the calculation including NN final 
state interaction (dashed curves) and the calculation with additional irN rescattering 
(solid curves). The left, middle and right panels represent the total cross section for 
jd — > ir~pp, ir + nn and ir°np, respectively. 



Fig. 14.11 shows that the differences between the full and the spectator model calcu- 
lations clearly demonstrate the importance of rescattering effects, in particular for 
the 7T° channel. One can see that final state interactions lead to a strong reduction 
of the total cross section. In the case of charged pion photoproduction reactions, the 
final state interaction effects are small in comparison with the case of neutral pion 
photoproduction reaction. The differene between the solid and dashed curves in Fig. 
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14. II shows that the nN rescattering changes the final results only by a few percent. A 
possible explanation for this smallness comes from the fact that the scattering length 
of the most important 5*-wave 7riV-scattering is about two orders smaller than the 
one of iViV-scattering. In the energy range of the A(1232) resonance, one finds the 
strongest manifestation of rescattering effects. In the case of tt~ and 7r + channels, 
one sees that the NN and irN rescattering effects reduce the total cross sections for 
both channels in the energy range of the A(1232) resonance by about 5%. For lower 
and higher energies, the influence of rescattering effects changes the results for total 
cross sections of the charged pion channels with a few percent. With respect to 7r° 
channel, we see from the right panel of Fig. 14.11 that, the influence of rescattering 
effects is much more significant than the case of charged pion channels. It leads to a 
strong reduction of the total cross section by about 35% in the energy range of the 
A(1232) resonance. 




Figure 4.2: The ratio <J tot as a function of the photon energy in the labora- 

tory frame. The left, middle and right panels represent the ratio for 'yd — > TT~pp, ir + nn 
and ir°np, respectively. 



In order to show in greater detail the relative influence of rescattering effects on the 
total cross sections, we show in Figs. 14.21 and 14.31 the ratios cr^ +NN+nN / a\^ t and 
<J tot +NN+nN l a tot +NN \ respectively, as a function of the photon energy in the labora- 
tory frame. Here denotes the total cross section in the impulse approximation, 
a toi +NN denotes the total cross section including NN rescattering and of ^ +Ar7V+7rAr 
denotes the total cross section with additional nN rescattering (see section l3~3*j) . 

Fig. 14.21 shows the combined effect of NN and nN rescattering as a function of 
the photon energy. This emphasize the important role of final state interactions 
for the total cross sections, especially for tt° photoproduction on the deuteron. It 
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Figure 4.3: The ratio <r^ +NN+nN fer^ +NN as a function of the photon energy. The left, 
middle and right panels represent the ratio for — > vr~pp, ir + nn and Tr°np, respectively. 



is obvious that for charged pion channels the ratio has a similar shape, but in the 
case of 7T° production, a totally different shape is observed. At energies near to 
the pion threshold, the effect of rescattering increases the total cross section for 7r + 
production than for 7r~ production. For tt + production we found that at 150 MeV 
the ratio is 1.32, while it reads 1.01 for tt~ production at the same energy. For tt° 
production, the situation is completely different since final state interactions lead 
to a strong reduction which amounts to about 60% at 200 MeV. This demonstrates 
that rescattering effects play a very important role, in particular, in the case of 7T° 
photoproduction. 

As next we show the relative effect of ttN rescattering in Fig. 14.31 In general, we 
obtain a similar shape as in the case of ff^ 4 "™" 1 "^/^ (see Fig. l4.2|) . especially that 
the peaks have approximately the same energy. Fig. 14.31 shows that the influence 
of ttN rescattering in the region of the A resonance is much less important than 
NN rescattering. For energies near to the pion threshold, we see that the effect of 
ttN rescattering increases the total cross sections for charged pion channels. This 
increment is noticeable in the case of ir + production since the ratio reads 1.033 at 
150 MeV. Furthermore, we see that the effect of ttN rescattering in the case of 
7T° production leads to a strong reduction for the total cross section, in particular 
at 200 MeV. A very big effect of the ttN rescattering in the threshold region was 
also found in Refs. [HJ E01 E2] • As already mentioned in that the charged 
pion channels are also of importance in the threshold region, because a big n^N 
rescattering effect is certainly possible. 
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4.1.1 Comparison with Experimental Data 

Here we compare our results for the total cross sections of the reaction 'yd — > nNN 
including NN and nN rescattering in the final state with experimental data. In 
our comparison with experiment we concentrate our discussion on 7r~ and tt° photo- 
production on the deuteron, since data for ir + production in the A(1232) resonance 
region are not available. In the case of ix~ production we compare our results with 
the experimental data from Refs. [2B E3] ■ For tt° production we compare our 
results with the experimental data from |22j . 
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Figure 4.4: Total cross sections for tt~ (left panel) and 7r° (right panel) photoproduction 
on the deuteron. Solid curves show our results with NN and ttN rescattering in the final 
state and dashed curves show the results in the impulse approximation. The experimental 
data are from [2J (ABHHM), [S3] (Frascati) and [SI] (Asai) in case of n~ production and 
from [22] (TAPS) in case of 7r° production. 



Fig. 14.41 shows our results for the total cross sections for 7r~ and 7r° photoproduc- 
tion on the deuteron compared with experimental data. One readily notes, that the 
spectator approach can not describe the experimental data in the A(1232) resonace 
region (see also [23]), especially in the case of n° production. Note, that the con- 
tribution of final state interactions in the case of 7r° production is much larger than 
that for 7T" production. Furthermore, final state interaction effects lead to a strong 
reduction of total cross section. The inclusion of such effects improve the agreement 
between experimental data and theoretical predictions. Indeed, a quite significant 
contribution from NN rescattering is found. Only in the center of the peak for tt° 
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production our model overestimates the measured total cross section by about 6%. 
Moreover, we found that the ttN rescattering changes the final results only by a few 
percent. As mentioned before, this smallness can be explained by the fact that the 
most important S*-wave ttN scattering length is about two orders smaller than the 
one for NN scattering. 

These results clearly show that the rescattering effects are significant and reduce 
the total cross section in the A(1232) resonance region. This means, in particular 
with respect to a test of theoretical models for pion production amplitudes on the 
neutron, that one needs a reliable description of the rescattering process. Com- 
pared to experimental data, one readily finds that the sizable discrepancies without 
rescattering are largely reduced and that a reasonable agreement with the data is 
achieved. 

4.1.2 Comparison with other Theoretical Predictions 

Pion photoproduction on the deuteron in the impulse approximation has been stud- 
ied by Schmidt et al. [23J neglecting all kinds of final state interactions and other 
two-body operators. Since data for n + and n° photoproduction on the deuteron 
were absent at the time of these calculations, the authors had no possibility to com- 
pare their predictions with experimental data for these channels. A comparison with 
experimental data was possible only for ir~ production, where they found a slight 
overestimation of the data. They reported that the reason for that is an overesti- 
mation of the elementary reaction on the neutron. They mentioned also that the 
differences between their theoretical predictions and the experimental data show 
very clearly that the calculation of pion photoproduction on nuclei in the nucleon 
spectator model can only be considered as a first step towards a more realistic de- 
scription of this process. As noted in ^21 I2H 1213 E2] the effect of NN rescattering 
is important in the incoherent pion photoproduction on the deuteron, especially for 
small pion angles. 

Laget [121 EH 120] and Blomqvist and Laget [TT] have investigated pion photopro- 
duction on the deuteron with the inclusion of pion rescattering and NN final state 
interaction within a diagramatic ansatz. They used the elementary photoproduction 
operator of Blomqvist and Laget ^Tj as input in their calculations. At the time of 
these calculations, a comparison with experimental data was possible only for 7r~ 
production since data for ir + and 7r° production in the A(1232) resonance region 
were not available. The agreement of their predictions including final state interac- 
tions with the experimental data of the cross sections of the reaction jd — ► ir~pp 
is quite good. They found that the final state interaction effects are small for the 
charged pion photoproduction reactions in comparison with the case of the neutral 
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channel. As mentioned in [22] , these predictions are significantly above the data for 
7t° photoproduction on the deuteron. A possible reason for this may be that they 
used the Blomqvist and Laget parametrization ^1] of the elementary photoproduc- 
tion amplitude. This amplitude is constructed using different A parametrizations 
for neutral and charged pion photoproduction. It gives a satisfactory fit to the am- 
plitude for charged pion photoproduction, but it is not able to describe the neutral 
pion photoproduction from the proton. An attempt to remedy this defect in [55J led 
to a 7T° photoproduction amplitude which is not very suitable for the use in nuclear 
calculations. 

Levchuk et al. [24J studied quasifree tt° photoproduction from the neutron via the 
d(7, ir°)np reaction using the elementary photoproduction operator of Blomqvist 
and Laget [TT] . The contributions from the pole diagrams as well as one-loop dia- 
grams both with np and nN rescattering were taken into account. They wanted to 
explore the possibility of measuring the E\ + /Mi + ratio via photoproduction from 
the quasifree neutron. The isospin 1 = 3/2 component of this ratio characterizes 
the relative strength of the recently much discussed (see e.g. [2EI EZ|) quadrupole 
£'2-excitation of the A resonance. The idea was that the 71(7, 7r°)n reaction would be 
very useful for the isospin separation of the multipoles. In agreement with the results 
from the Laget model j!2j . Levchuk et al. find that the largest effects disturbing 
the extraction of the multipoles for quasifree neutrons arise from the np final state 
interaction. They predict that these effects lead to a strong reduction of the cross 
section at pion forward angles, but are much less important for backward angles. 
They found also that the correction due to np rescattering decreases with increasing 
pion angle and becomes to be less than 8% at 9 W > 90°. Furthermore, they pointed 
out that the contribution of the proton pole diagram and the one of ttN rescattering 
are negligible. The experimental data from [22] for the ^(7, 7r°)np reaction qualita- 
tively support this prediction since the disagreement with the spectator approach is 
most severe at pion forward angles but less pronounced at backward angles. How- 
ever, a comparison of the data to the Laget model including final state interactions 
shows some unexplained reduction of the cross section at backward angles. 

Recently and during the calculations of this thesis work, Levchuk et al. |2S] studied 
the inclusive reaction (£(7, ir)NN in the A resonance region. This calculation is based 
on the use of the diagramatic approach. Pole diagrams and one-loop diagrams with 
NN and ttN rescattering in the final state are considered. The authors in [23] 
pointed out that the main difference between their calculation and the one of Ref. 
[T2~] is that a more realistic version of the elementary pion photoproduction operator 
is used. It is taken in on-shell form and calculated using the SAID 28 s and MAID J7j 
multipole analyses. This operator is constructed in the c. m. frame. Therefore, 
it has to be transformed to an arbitrary frame of reference to be used as input 
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in calculations on light nuclei. This may be done by a Lorentz boost of all four 
momenta on which the elementary amplitude depends. Their predictions for total 
and differential cross sections including final state interactions show good agreement 
with the experimental data. 

Now, we compare our results for the total cross section of the processes (IH.lj) with 
the theoretical predictions from [23] as shown in Fig. 14.51 



300 



300 




800 



■a 



b 



150 250 350 450 
co^ b [MeV] 



150 250 350 450 
co^ b [MeV] 



"i r 



7t°np 




150 250 350 450 
(0y ab [MeV] 



Figure 4.5: Total cross sections for pion photoproduction on the deuteron obtained within 
the impulse approximation (dotted curves) in comparison with the calculation when NN 
and ttN final state interactions are taken into account (solid curves). The results of 
Levchuk et al. [2^] as shown by the dash-dotted curves for the impulse approximation and 
by the dashed curves with NN- and 7riV-rescattering. 



In the impulse approximation, Fig. 14.51 shows that our results for charged pion 
photoproduction (dotted curves) are higher than the results including rescattering 
effects (solid curves). This means that rescattering effects lead to a reduction of 
the total cross sections in our case. In the same figure, we see from the difference 
between the dash-dotted and dashed curves that the theoretical predictions from [25J 
display the opposite situation. In the next section we will give an explanation for this 
opposite situation in some detail. In the case of tt° photoproduction, rescattering 
effects lead in both cases to a strong reduction of the total cross section. 

Including final state interactions we found small discrepancies between our calcula- 
tions which are given by the solid curves in Fig. 14.51 and the predictions of Ref . {213 
which are given by the dashed curves in the same figure. A reason for that may 
be the use of different pion photoproduction operators. As mentioned in section 
13.3.21 that we used in our calculation the separable representation of the Paris NN 
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potential ^ElIM! to solve the Lippmann-Schwinger equation for the NN scattering. 
In the authors used another type of NN potentials which is the Bonn OBE 
potential model (OBEPR) 03]. This choice of different NN potential models may 
also be a reason for the small discrepancies between both theoretical predictions. 
In the case of 7rN scattering, the separable energy-dependent ttN potential built in 
[TB] is used in both theoretical predictions. 
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4.2 Differential Cross Section 

In this section we discuss our results for the differential cross sections and compare 
with experimental data and other theoretical predictions. We begin with presenting 
our results in Fig. 14.61 




Figure 4.6: Differential cross sections for pion photoproduction on the deuteron within the 
impulse approximation (dotted curves) in comparison with the calculation including NN 
final state interaction (dashed curves) and the calculation with additional irN rescattering 
(solid curves). The left, middle and right panels represent the differential cross section for 
ryd — > Tt~pp, ir + nn and 7r°np, respectively. 



As in the case of total cross section in the previous section, Fig. 14.61 shows that the 
difference between the full and the spectator model calculations clearly demonstrate 
the importance of rescattering effects. One sees also that ttN rescattering changes 
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the final results for the differential cross sections only by a few percent. Furthermore, 
for charged pion channels, one can see that the effect of final state interactions 
becomes small for low and high energies, but it becomes maximal at energies near 
to the A(1232) resonance. In the case of 7r° channel, it is obvious that the effect of 
rescattering is important for all energies, in particular for forward pion angles. In 
the case of both charged and neutral pions, we see that rescattering effect is quite 
small for backward pion angles. Looking at the right and middle panels in Fig. 14.61 
one can see that at 6 n = 0, the differential cross section increases with increasing 
photon energy. This increase comes mainly from the Born terms. We found that 
more than 70% from the values of the differential cross section at 9 n = comes 
from the Born terms and less than 30% comes from the contribution of the A(1232) 
resonance (see Fig. 14. 7|) . 
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Figure 4.7: Differential cross section for 7r~ photoproduction on the deuteron in the 
impulse approximation at different photon energies. The solid curve shows the results 
using only the Born terms and the dotted one shows the contribution of the A(1232) 
resonance. 



In order to show in more detail the relative importance of final state interac- 
tions on the differential cross sections, we show in Figs. 14.81 and 14.91 the ratios 

da IA+NN+irN j d(J IA and d(J IA+NN+*N J folA+NN ? respective l y; as a f unctl ion of the 

pion angle in the laboratory frame. 

Fig. 14.81 shows the relative effect of both NN and nN final state interactions by 
the ratio da IA+NN+nN /da IA as a function of the pion angle. As mentioned in the 
previous section, we see that the ratio has a similar shape for charged pion channels. 
It is also obvious that rescattering effects are important for forward pions and much 
less important for backward pions. In the case of 7r° production, we see that the 
influence of final state interaction effects on the differential cross section is significant. 
Rescattering effects lead to a strong reduction of the differential cross sections, in 
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particular for small pion angles. The contribution decreases quickly when the pion 
angle increases and becomes very small for backward pion angles. 
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Figure 4.8: The ratio da IA+NN+nN /da IA as a function of the pion angle in the laboratory 
frame at different photon energies. The left, middle and right panels represent the ratio 
for the jd — > ir~pp, ir + nn and 7r°np, respectively. 



Now, we would like to clarify the relative effect of irN rescattering in Fig. 14.91 bv the 
ratio da IA+NN+nN /da IA+NN . Here one can also see that we obtain a similar shape 
for the charged pion channels, but a totally different shape for the 7r° channel is seen. 
The effect of irN rescattering is seen much less important than NN rescattering. 
This effect is less than 1% for charged pion channels and less than 5% for neutral 
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channel. In the energy range of the A(1232) resonance, we find that this effect is 
maximal at pion forward angles. 
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Figure 4.9: The ratio da IA+NN+7rN /da IA+NN as a function of the pion angle at different 
photon energies. The left, middle and right panels represent the ratio for the jd — * ir~pp, 
ir + nn and 7T°np, respectively. 



4.2.1 Comparison with Experimental Data 

Here we compare our results for the differential cross sections with the experimental 
data from [21] for 7r~ production and from [22] for 7T° production. This comparison 
is shown in Fig. 14.101 As mentioned in the previous section, there are no data 
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available for 7r + production on the deuteron in the A(1232) region so that we can 
not compare our predictions with experimental data for this channel. Therefore, 
we concentrate our discussion on ir~ and 7r° production. The effect of final state 
interactions in the case of charged pion channels is expected to be quite different in 
comparison with the neutral channel. 
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Figure 4.10: Differential cross sections for 7r~ (left panels) and 7r° (right panels) photopro- 
duction on the deuteron. Solid curves show our results include NN and ttN rescattering 
and dashed curves show the results in the impulse approximation. The experimental data 
are from |21| for tt~ production and from \2'2\ for 7T° production. 

In Ref. [22] the differential cross sections for the reaction di^y, n )np are given in 
the so-called jN c. m. frame. Therefore, in order to compare our results for the 
differential cross sections of the reaction <i(7, 7i°)np with the experimental data from 
[2*2"] we need to transform the differential cross sections from the rest frame of the 
deuteron to the jN c. m. frame. The Jacoby determinant which we need for this 
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purpose is given in appendix El The pion angle in the 'yN c. m. frame is denoted 
by 6^ 

Fig. 14.101 shows that the spectator approach can not describe the experimental data 
for differential cross sections. One sees also that the effect of final state interactions 
leads to a strong reduction of the differential cross sections, especially for forward 
pion angles but it is much less important for backward angles. The inclusion of such 
effects improves the agreement between experimental data and our calculations. 
After including final state interactions we obtain in general a good agreement with 
the experimental data for differential cross sections. Small discrepancies are found 
only at backward pion angles. We see that the agreement with the experimental 
data of the reaction — > 7r~pp measured in a bubble chamber experiment |2~T] . is 
quite good. The comparison with the TAPS data j22j for the reaction jd — > ir°np 
yields a quite reasonable description of the experimental differential cross section, 
in particular in the energy range of the A(1232) resonance. The effect of final state 
interactions is found to be smaller for charged pion channels than in the case of 
neutral channel. 
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Figure 4.11: Differential cross section for ir° photoproduction on the deuteron as a func- 
tion of pion angle in the c. m. frame. Solid curves show our results include NN and 
ttN rescattering and dashed curves show the results in the impulse approximation. The 
experimental data are from |22j . 

As shown in Fig. 14.41 that our predictions for the total cross section overestimate 
the experimental data for 7r° photoproduction in the energy region between 340-360 
MeV by about 6%. In order to see how good is the agreement with the experimental 
data for the differential cross section at this energy region, we predict in Fig. 14.111 
a comparison between our results and the experimental data from TAPS |22j . It 
is noticeable that our predictions overestimate the experimental differential cross 
section at this energy, is particular at backward pion angles. 
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4.2.2 Comparison with other Theoretical Predictions 

In this section we compare our results for the differential cross sections with the 
theoretical predictions of Levchuk et al. [25J. As in [2H], the pion angle is given in 
the 7iV c. m. frame in the case of 7r° channel. 

It is worthwhile first to point out that, in the impulse approximation Fig. 14.121 
shows the following interesting features. Without final state interactions we found 
a big difference between our calculation and the prediction from [2EJ for charged 
pion channels at pion forward angles (see the difference between the dotted and the 
dash-dotted curves in Fig. 14.12)1 . Therefore, we checked in this work where this big 
difference comes from. 

First, we assume that the authors in Ref. [23] used a wrong NN antisymmetrization. 
That is why we examine using a wrong NN antisymmetrization procedure in our 
work. We found that if we use a wrong NN antisymmetrization for the s = 
channel, i.e., 

V2> (1 ¥i> (2 0l s ™*>l*™*>> 



V2} {1) \Pi) {2) )\sm s }\tm t }, (4.2) 

we obtain results as shown in Fig. I4.1HI (dotted curves) with the same energy shape 
as that of Levchuk et al. j2S] (dashed curves). This probably means that, they did 
not use the correct NN antisymmetrization in their calculations. 

Moreover, since we use a coupled basis state for the NN antisymmetrization in our 
calculations, we also checked the use of an uncoupled, i.e., helicity basis which was 
used in We found that both approaches give the same results. This confirms the 
NN antisymmetrization which we use and one can only suspect that the difference 
to the results of Ref. |2S] may originate from an error in the antisymmetrization. 

To investigate the influence of 7riV — > KY rescattering in kaon photoproduction on 
the deuteron, the authors in (SU used the same procedures as in Levchuk et al. j2Sj 
to calculate pion photoproduction on the deuteron in the impulse approximation. 
They obtained results as shown in Fig. 14.141 (dashed curves). It is obvious that 
their results using the elementary production operator from MAID [17 analysis are 
in good agreement with our results (solid curves) using the elementary production 
operator presented in chapter El It is also clear that in both cases the differential 
cross section has a high value at 6 n = and not as in the case of Ref. |2Sj- This 
suspects that the model in Ref. has a kind of error. 



\pip 2 ,sm s ,tm t ) = (|pi> (1) |P2> (2) + (-) 
and a correct one for s = 1 channel, i.e., 
\pip 2 ,sm s ,tm t ) = —= (\pi) {1) \p2) {2) ~ (-) 
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Figure 4.12: Differential cross sections for pion photoproduction on the deuteron in com- 
parison with the results from |25| at different photon energies. Full curves show our 
results when rescattering effects are included and dotted curves show our results in the 
impulse approximation. Dashed and dash-dotted curves show the results from with 
and without rescattering effects, respectively. 



After including final state interaction effects, Fig. 14.121 shows that the agreement of 
our results (solid curves) and the theoretical predictions from [2E] (dashed curves) 
is good. A small discrepancy between both calculations is observed. A reason for 
that may be due to the use of different pion photoproduction operators. 
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Figure 4.13: Differential cross section for ir~ photoproduction on the deuteron in the 
impulse approximation. The full curve shows the results of our calculation and the dashed 
one shows the results from |25| . The dotted curve shows our results with a wrong NN 
antisymmetrization for s = channel, as given in Eq. (|4.1jl . 
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Figure 4.14: Differential cross section for ir~ photoproduction on the deuteron in the 
impulse approximation. The full curve shows our results using the elementary production 
operator of chapter [21 and the dashed one shows the results from |56| using the elementary 
reaction amplitude from MAID [EI analysis. 
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4.3 Spin Asymmetry and Gerasimov-Drell-Hearn 
Sum Rule 

Recently, several experiments to test the GDH sum rule both on the proton and 
on the neutron are now in preparation or planned at different laboratories around 
the world (MAMI, ELSA, LEGS, GRAAL and TJNAF). This makes the study of 
polarization observables nowadays of great interest in the field of intermediate energy 
nuclear physics. The spin asymmetry of the total photoabsorption cross section, 
entering into the Gerasimov-Drell-Hearn (GDH) sum rule |H7| I58j. is of particular 
interest. Therefore, we evaluate in this section the contribution of incoherent single 
pion photoproduction including NN and nN rescattering to the spin asymmetry 
and the corresponding GDH sum rule for the deuteron using the model developed 
in chapter El 

The verification of the GDH sum rule is an important issue, which enables a check of 
some fundamental physical principles related to the spin of the nucleon. The GDH 
sum rule connects the anomalous magnetic moment of a particle with the energy 
weighted integral from threshold up to infinity over the spin asymmetry of the 
total photoabsorption cross section, i.e., the difference of the total photoabsorption 
cross sections for circularly polarized photons on a target with spin parallel and 
antiparallel to the spin of the photon (see Fig. 14.15]) . In detail it reads for a particle 
of mass M, charge eQ, anomalous magnetic moment k and spin S 




Figure 4.15: Illustration of the relative spin orientation of the incoming photon and the 
target nuclei in the GDH sum rule. 

I GDH M = = j~ d ±{a p (k)-a\k)) , (4.3) 

where a p ^ A (k) denote the total absorption cross sections for circularly polarized 
photons on a target with spin parallel and antiparallel to the photon spin, respec- 
tively, and the anomalous magnetic moment is defined by the total magnetic moment 
operator of the particle 



M = (Q + k)4-S. 



(4.4) 
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This sum rule gives a very interesting relation between a magnetic ground-state 
property of a particle and an integral property of its whole excitation spectrum. In 
other words, this sum rule shows that the existence of a nonvanishing anomalous 
magnetic moment points directly to an internal dynamic structure of the particle. 
Furthermore, because the left-hand side of Eq. ()4.3|) is positive, it tells us that the 
integrated, energy-weighted total absorption cross section of a circularly polarized 
photon on a particle with its spin parallel to the photon spin is bigger than the 
one on a target with its spin antiparallel, if the particle posesses a nonvanishing 
anomalous magnetic moment. 

Since proton and neutron have large anomalous magnetic moments (k p = 1.79 and 
K n = —1.91), one finds correspondingly large GDH sum rule predictions for them, 
i.e., I^ DH (oo) = 204.8 /ib for the proton and I^ DH (oo) = 233.2 /ib for the neutron. 
Applying the GDH sum rule to the deuteron, one finds a very interesting feature. 
The deuteron has isospin zero, ruling out the contribution of the large nucleon 
isovector anomalous magnetic moments to its magnetic moment. Therefore, one 
expects a very small anomalous magnetic moment for the deuteron. In fact, the 
experimental value is = —0.143 resulting in a GDH prediction of lf DH (oo) = 
0.65 /A), which is more than two orders of magnitude smaller than the nucleon values. 

In recent years, Arenhovel et al. [59J have evaluated explicitly the GDH sum rule 
for the deuteron by integrating the difference of the two total photoabsorption cross 
sections with photon and deuteron spins parallel and antiparallel up to a photon 
energy of 550 MeV. There, three contributions are included: (i) the photodisinte- 
gration channel 'jd — > np, (ii) the coherent pion photoproduction 7<i — > TT°d and (iii) 
the incoherent single pion photoproduction — *> ttNN '. In their calculation of the 
7<i — > tcNN contributions to the GDH integral, the authors restricted themselves 
to the impulse approximation using the spectator nucleon approach of Schmidt et 
al. [23J. For the total GDH value from explicit integration up to 550 MeV, they 
found a negative value If DH (550 MeV) = — 183/ib. However, some uncertainty lies 
in the contribution of the incoherent single pion production channel because final 
state interactions and other two-body operators are not included in the spectator 
nucleon model of Ref. [23J. As discussed previously, final state interactions in inco- 
herent single pion photoproduction on the deuteron play an important role on the 
differential and total cross sections. Therefore, the influence of rescattering effects 
on the spin asymmetry and the corresponding GDH sum rule has to be investigated. 

In this work, the influence of NN and nN rescattering effects in incoherent single 
pion photoproduction to the spin asymmetry and the GDH sum rule for the deuteron 
is investigated. We evaluated explicitly the GDH sum rule for the deuteron by 
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integrating the difference of the two total photoabsorption cross sections with photon 
and deuteron spins parallel and antiparallel up to a photon energy of 550 MeV. The 
upper integration limit of 550 MeV is chosen, because we consider only single pion 
production. The contributions which arise from coherent pion photoproduction and 
photodisintegration channels are given in Ref. [59J to which we refer for more details. 

Our results for the spin asymmetry are collected in Fig. 14.161 for the individual 
contributions from the different charge states of the pion. Fig. 14.161 shows our 
results for the total photoabsorption cross sections for circularly polarized photons 
on a target with spin parallel a p (upper part) and antiparallel a A (middle part) to 
the photon spin. In the lower part of Fig. 14.161 the difference of the cross sections 
/ -a A for incoherent pion photoproduction on the deuteron including final state 
interactions is shown. For comparison, we also show our results on the free nucleon 
by the dash-dotted curves in Fig. 14.161 using the elementary production amplitude 
constructed in chapter 121 One notes qualitatively a similar behaviour for the spin 
asymmetry although the maxima and minima are smaller and also slightly shifted 
towards higher energies for the deuteron. The bottom panels in Fig. 14.161 show also 
that final state interactions lead to a strong reduction of the spin asymmetry in the 
energy region of the A(1232) resonance. This reduction becomes more than 35% 
for 7T° photoproduction and becomes more than 15% for charged pion channels. It 
is also obvious that a p is much greater than a A , in particular in the case of n° 
production. 

The results for the GDH sum rule are depicted in Fig. 14.171 and Table 14.11 for the 
individual contributions from the different charge states of the pion. Their total 
sum to the GDH sum rule is shown in Fig. 14.181 Our results on the free nucleon are 
also shown for comparison (see the dash-dotted curves in Figs. 14.171 and I4.18j) . It 
is obvious that a large positive contribution to the GDH sum rule comes from the 
7T° production channel whereas the charged pions give a negative but - in absolute 
size - smaller contribution to the GDH value. Up to an energy of 550 MeV one 
finds for the total contribution of the incoherent pion production channels a value 
/ 7 G £W550MeV)=87/ib. 

A very interesting and important result is the large negative contribution from the 
7r ± channels and the large positive contribution comes from the 7r° channel to the 
GDH value. Hopefully, this low energy feature of the GDH sum rule could be checked 
experimentally in the near future. At the same time, precise data on a p — a A from 
a direct measurement is urgently needed. 

Last but not least, we would like to conclude that the results presented here for the 
spin asymmetry and the GDH sum rule can be used as a basis for the simulation 
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Table 4.1: Contributions of incoherent single pion photoproduction to the GDH integral 
for the deuteron integrated up to 550 MeV in /ib. 



rpartinn T GDH T GDH T GDH 

reaction i IA J-ia+nn 1 ia+n 



-fd^ppiT- -73 -87 -88 
7 d -> nnvr + -27 -39 -41 
^d-^np-K 287 220 216 



-fd-^nNN 187 94 87 



of the behaviour of polarization observables and for an optimal planning of new 
polarization experiments of the reaction jd — > nNN. 
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Figure 4.16: Total absorption cross sections for circularly polarized photons on a target 
with spin parallel a p (upper part) and antiparallel a A (middle part) to the photon spin. 
Lower part shows the difference of the cross sections. Dotted curves show the results 
in the impulse approximation, dashed curves show the results including AW rescattering 
and solid curves show the results with additional ttN rescattering. The dash-dotted curves 
show the results for ir~ on the neutron (left panels), tt + on the proton (middle panels) 
and 7T° on both the proton and the neutron (right pannels). 
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Figure 4.17: The Gerasimov-Drell-Hearn integral (see Eq. (|4.3jl ) as a function of the upper 
limit of integration for the different channels of incoherent single pion photoproduction on 
the deuteron and the nucleon. Notation of the curves as in Fig. 14.161 




Figure 4.18: Summation of the contributions of the three channels of incoherent single 
pion photoproduction to the GDH sum rule for the deuteron and the nucleon as function 
of the upper integration energy. See caption of Fig. 14.161 for meaning of the curves. 



Chapter 5 

Summary and Outlook 



5.1 Summary 

The subject of this work was the investigation of the influence of final state in- 
teraction effects in incoherent single pion photoproduction on the deuteron in the 
A(1232) resonance region. The transition matrix elements are calculated in the 
frame of time-ordered perturbation theory, using an elementary production opera- 
tor, which yields a good agreement with experimental data, and including final state 
interaction effects. 

For the elementary 7iV — > ttN photoproduction amplitude we used the effective La- 
grangian model of Schmidt et al. [23J. This amplitude contains besides the standard 
pseudovector Born terms the resonance contribution from the A(1232) excitation. 
The Born terms of the elementary production amplitude are determined in pseu- 
dovector tcN coupling and supplied with a form factor. The A(1232) resonance is 
considered in both the s- and w-channel. The parameters of the A resonance and the 
cut-off of the form factors are fixed on the leading photoproduction multipole am- 
plitudes. We found a good agreement between our results for total and differential 
cross sections and experimental data. 

This elementary production operator is then used to study pion photoproduction 
on the deuteron. In addition to the impulse approximation, where all kind of final 
state interactions and other two-body operators are neglected, we have included 
as presumably dominant final state interaction effects two-body interactions in the 
NN- and 7riV-subsystems. As models for the interaction of the NN- and ttN- 
subsystems we used separable interactions which are fitted to the phase shift data for 
NN and ttN scattering. For NN scattering we have included all partial waves with 
total angular momentum J < 3 by numerical solution of the Lippmann-Schwinger 
equation [17| with the separable representation of the Paris NN potential jlHl HH] ■ 



67 



68 



Chapter 5. Summary and Outlook 



In the case oiirN rescattering we have considered the S-, P- and D- waves by solving 
the Lippmann-Schwinger equation for a separable energy- dependent nN potential 
from ITS]. 

For the total cross section we found that the influence of NN and ttN rescattering 
effects on the total cross section is significant. These effects reduce the total cross 
sections for the charged pion photoproduction reactions in the A(1232) resonance 
region by about 5%. In the case of 7T° photoproduction reaction, rescattering is much 
more important, reducing the total cross section by about 35% on the maximum. 
Furthermore, it was found that ttN rescattering is much less important compared 
to iViV rescattering. In comparison with experimental data, we found that the in- 
clusion of NN and irN rescattering effects leads to an improved agreement between 
experimental data and theoretical predictions. Only in the maximum of 7r° pro- 
duction our model overestimates the measured total cross section by about 6%. In 
comparison with the theoretical predictions from we found that for charged 
pion production the inclusion of rescattering effects leads to a reduction of the total 
cross sections in contrast to the theoretical predictions from [25J. This different be- 
haviour may have its origin in the difference already for the impulse approximation 
for which [23] predicts a 14% lower total cross section. In the case of 7r° channel, 
also in [25j a strong reduction by rescattering has been found. After including final 
state interactions we found, however, only small differences between our calculations 
and the predictions from |25j . 

In the case of the differential cross section one may draw the following conclusions. 
The inclusion of NN and irN final state interactions reduces the differential cross 
sections for the charged pion channels mainly at pion forward angles by about 15%. 
The reduction is much stronger for the 7r° channel by about 40% at pion forward 
angles. At pion backward angles the influence is much less important. As already 
noted for the total cross section, the irN rescattering changes the final results only 
by a few percent. In comparison with experiment, after including final state in- 
teractions we have obtained a satisfactory agreement with the experimental data. 
Small discrepancies were found only at backward pion angles. In comparison with 
the theoretical predictions from |2Sj, without final state interactions we found a big 
difference between both calculations for the differential cross sections for charged 
pion channels at pion forward angles. After including rescattering effects we have 
obtained a good agreement with the theoretical predictions from |25j . 

Finally, we have evaluated in this work the contribution of incoherent single pion 
photoproduction to the spin asymmetry of the total photoabsorption cross section. 
We found that the inclusion of final state interactions leads to a strong reduction of 
the spin asymmetry in the energy region of the A (1232) resonance. This reduction 
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amounts to about 35% for 7r° photoproduction and about 15% for charged pion 
channels. The corresponding GDH integral was also evaluated up to 550 MeV. We 
found that a large positive contribution to the GDH sum rule comes from the ir° 
production channel whereas the charged pions give a negative but - in absolute size 
- smaller contribution to the GDH value. Up to an energy of 550 MeV we have 
obtained for the total contribution of the incoherent pion production channels a 
value 7^^(550 MeV) =87 /xb. 

5.2 Future Extensions 

The studies we have discussed here will serve as the basis for further investigations 
including the dynamics of the ttNN system in a more complete way. Undoubtedly, 
there is still a lot of work to be done both experimentally and theoretically. In- 
deed, many challenging and interesting lessons have yet to be learned before a deep 
understanding of the incoherent pion photoproduction process on the deuteron will 
emerge. In first instance, this work could be continued by the further refinement of 
the elementary production operator. Modifying the elementary pion photoproduc- 
tion operator on the free nucleon above the two pion threshold with, for example, 
the inclusion of the u and p meson exchange in the t-channel and contributions 
at higher energies may improve our results for the spin asymmetry and the corre- 
sponding GDH sum rule for the deuteron. Our goal is to have an operator that can 
describe the elementary process on the free nucleon reasonably well over a larger 
energy region and that is suitable for applications to nuclear systems. 

In view of the importance of the first order rescattering it is natural to ask about 
the role of higher order rescattering terms. Answering this question points to the 
necessity of a three-body approach, where the final state interaction is included to 
all orders. This may result in a better agreement between experimental data and 
theoretical predictions for 7r° photoproduction. Furthermore, the inclusion of such 
rescattering terms may be important in studying polarization observables. Studying 
these observables will give us much more detailed information and thus will provide 
much more stringent tests for theoretical models. In fact we plan to embark on such 
a study in the near future. 

In the long run, one would need also to extend the formalism to the threshold region 
for which the elementary production operator has to be improved. This process is 
of great interest since experimental data for the reaction d^j, 7r°n)p have been mea- 
sured recently in Mainz (MAMI/TAPS) and Saskatoon (SAL) [60J. Moreover, the 
formalism should be extended to investigate coherent and incoherent elect reproduc- 
tion of pions on the deuteron including final state interaction effects in both the 



70 



Chapter 5. Summary and Outlook 



threshold and the A(1232) resonance regions in order to analyze recent results from 

MAMI inn. 
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Appendix A 

General Notations and 
Conventions 



A.l Dirac Algebra 



Throughout this work we adopt the natural system of units where K — c = 1; 
conversion factor: He = 197.32696 MeV.fm 1 . We also follow the convention of 
Bjorken and Drell 43J, where the contravariant space-time four- vector is defined as 



(x° , x 1 , x 2 , x 3 ) 



(t,x) = (t,x,y,z), 
and the covariant space-time four-vector is given by 
Xp = (x ,x 1 ,x 2 ,x 3 ) = (t,-x) = (t,-x,-y,-z) 

with the transformation matrix 
/ 1 \ 



0-100 
0-10 
\ -1 / 



(A.l) 



(A.2) 



(A.3) 



Likewise the contravariant four-momentum is 

V » EE (p°,p\p\p 3 ) = (E p ,p), 

and the scalar product between two four-momenta is given by 
p ■ q ee p^q^ ee E p E q - p-q. 



(A.4) 



(A.5) 



the calculations of this work we use units of MeV for energies, masses and momenta. 
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The Dirac matrices are 



( 7 °, 7), (A.6) 



with the matrix representation 

/ = r l ° ) ,? = ( ° f ) , (A.7) 



v o -ay" y — <t o 

where 70 is hermitian and 71, 72 and 73 are anti-hermitian. The Pauli matrices a = 
(a 1 , a 2 , a 3 ) are denoted by 

These matrices satisfy the anticommutation relations 

{a\a j } = <tV + ctV = 2%, (A.9) 
as well as the commutation relations 

[a\a j ] = crV-crV = 2ie ijk a k , (A.10) 
where represents the antisymmetric Levi-Civita tensor in R 3 . 

The Dirac matrices 7 satisfy the anticommutation relations 

= 7Y+7Y = 2<r. (ah) 

Important combinations of 7 matrices are the antisymmetric combination 

^ = ^[7 M ,7l (A.12) 



with components 

o K I " " U ! 



a* = ( °1 ° 1 and <7 W = i ( ° ) . (A.13) 



Other useful combinations are 

7 = 2 7 777 = 75 = ■fiie^peYl 11 = I j q 

ie-Wpal* = l5(-lulpla + 9upla + 9palu~ 9ualp) , (A.15) 

75^ = ^e^<v , (A.16) 
757. = -7.7s = ^W^V- (A.17) 
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The antisymmetric Levi-Civita tensor is defined by 

+1 for an even permutation (e.g. 0, 1, 2, 3) 



-fiupa 



-1 for an odd permutation 
if two or more indices are the same 



(A.18) 



The scalar product between 7 matrices and a four-momentum is written as 

7% = 7V-7-P = t- (A.19) 



The positive-energy four-component free Dirac spinor has the form 

/ Xs \ 

—* —* 

a ■ p 



u(p, s) 



E p + m\\ 



2m 



-Xs 



(A.20) 



and the negative-energy four-component Dirac spinor (antiparticle spinor) has the 
form 

v(p, s) 



/E p + m\ I 
V 2m ) 



a ■ p 



E p + m 

Xs 



Xs 



where E p = ^Jm? + p 2 . They are normalized as 

u(p, s)u(p, s) = 1 , 
v(p,s)v(p,s) = -1, 

where Xs are the two-component Pauli spinors with x+i 
1 , and the Dirac adjoint spinors are defined as 



u(p, s) 
v(p, s) 



(A.21) 



(A.22) 
(A.23) 



and X- 



(A.24) 
(A.25) 



Using Dirac spinors u and v, the Dirac equations may be written as 

{j> — m)u(p, s) = , 
{j) + m)v{p, s) = , 
which in terms of the adjoint spinors become 

u(p, s)($ - m) = 0, 
v(j), s)(J) + m) = 0. 



(A.26) 
(A.27) 

(A.28) 
(A.29) 
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A. 2 Normalization of States 

The general form of the state which describes the single particle X e {N, A,7r, 7} 
with momentum p and quantum numbers 2 a is \X;pa). It is also possible to use the 
form \p; a) when we know the type of the used particle. In case of more than one 
particle, the state has the form |X 1 ,X 2 ;pia 1 ,p 2 a 2)- 

For all the involved particles a covariantly normalized state is used in this work. 
The fermions are normalized as 

(N;fi\N;fi) = W 3 5 3 (P2 -Pi)^- (A.30) 

M N 

The normalized states for bosons are given by 



{n-q2\n-,qi) = (2n) 3 6 3 (q 2 -qi)2u; g , u q = y/m\ + f (A.31) 

{d-p 2 \d-p 1 ) = (2it)H\p-> 2 -n)2E d , E d =^Mj + fP (A.32) 

(7; k 2 \ r , ki) = {2n) 3 5 3 {k 2 - k^uj^ , uj y — \k\ — k . (A.33) 
The completeness relation for fermions then reads 

1 ^ / (A ' 34) 

and for bosons reads 

l =Sw?^ m - (A ' 35) 

A. 3 Spherical Basis 

— * 

For the construction of the physical pions from the pion field $ and for the evaluation 
of the spin and isospin operators, a and r, their representation in spherical basis is 
used. The basis vectors are given by 

= --^= (e x + i/jte y ) , for /j, = ±1 

e = e z . (A. 36) 

The spherical components of a given vector are determined by writing the scalar 
product 

a il = a-e il . (A. 37) 



2 These quantum numbers indicate about spin, isospin, spin projection etc. 
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If the spherical components of two vectors are given, then the scalar product of these 
two vectors in spherical basis is given by 

a-b = J2^(-Tb^ (A.38) 

and the vector product is given by 

(a x % = -iV2j2(^a lM6|l/x) • ( A - 39 ) 

Ma Mi) 
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Appendix B 



Useful Formulas for the 
Elementary Process 



In this appendix we give the useful formulas for the pion photoproduction reaction 
on the free nucleon. This process is schematically sketched in Fig. 12.11 Many of the 
formulas given here have been published elsewhere; our aim has been to collect all 
of the kinematic equations needed in this work. 

B.l Kinematics and Relevant Formulas 

The kinematics of pion photoproduction reaction on the free nucleon, ^yN — > ttN, are 
characterized by the four-momentum vectors k = (ou 7 , k) for the incident photon, 
Pi = (Ei,px) for the initial nucleon, p 2 = (E 2 ,P2) for the final nucleon and q = 
q) for the produced pion. Then the kinematical equations (energy-momentum 
conservations) are given by 

pi + k = p 2 + q, 




Ml 



(B.l) 



where Mn is the nucleon mass and m n is the pion mass. 



The usual Mandelstam kinematical variables are given by 



s = (P2 + q? = (Pi + k) 2 , 
t = (k-qf = (p 2 - Pl ) 2 , 
u = (pj, - kf = (p! - q) 2 , 



(B.2) 
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and they satisfy 
s + t + u = 



2M 2 N + ml 



(B-3) 




p 1= -k 



Figure B.l: Kinematics in the ttN cm. frame for the reaction 7JV — > irN. 



In the ttN cm. frame which given graphically in Fig. IB.II let W be the total energy 
of the system. Then 

s = W\ 

t = 2k ■ q — 2cu^ujg-+ m\ , 

u = -2k ■ q- 2uj^E 2 + M 2 N , (B.4) 
where p\ = —k and p 2 = —q in the cm. frame. The scattering angle 9 is given by 
(k-q) 



cos 9 



X . 



(B.5) 



(1*11*1) 

In the following, we will use the following notation to denote the 3-momentum in 
the system 

k 



k = \k\ k = -; q = I q \, q = - 
k q 



(B.6) 



The invariant mass of the ttN subsystem is given using the relativistic energy- 
momentum relations by 

2 

W 2 



M 2 N +p 2 + ^ml + q^ -(p 2 + q)' 



M 2 N + q 2 + ^ml + q 2 



M 2 N + 2M N J ab 



(B.7) 
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The various energies and momenta in the rrN cm. frame are given in terms of the 
invariant mass of the nN subsystem, W, by 

W 2 - Ml 

* = = <*,, (B.8) 



2W 



(W 2 -M 2 N + ml) 2 - Am 2 n W 2 



(B.9) 



W 2 -M 2 N + ml 



2W 



(B.10) 



Ei 



W 2 + M 2 



N 



2W 



(B.ll) 



Eo 



W 2 + M 2 N - ml 



2W 



[B.12) 



B.2 Multipole Decomposition of Amplitudes 

The on-shell X/j- matrix element is expanded in the cm. frame to four independent 
amplitudes as follows [3 



<g> | T fi (W) | fcg) 



AttW 



i(a ■ e)F 1 + (a ■ q)a ■ {k x e)F 2 + i(a ■ k)(q- e)F 3 

(B.13) 



M N L 

+i{B^q)(q-e)F A 



where Fi = Fi(W,x) with x = q- k are the well known CGLN-amplitudes. k and 
q are momentum unit vectors of the photon and meson, respectively. The isospin 
decomposition of the CGLN-amplitudes has the form 



p(-) [ r j> r o\ F (+) ( r M ' r o} , p(0) + 



(i = 1,2,3,4) 



2 2 

with the following commutator and the anti-commutator relations 

[r+,7b] = 2/ir; = (-r2/xr_ M) 
{r+,r } = (-) M {r-„r } = 25 M ,o , 



(B.14) 



(B.15) 
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where r M and T3 are the nucleon isospin matrices which act on the isospinor of the 
nucleon. 

With this convention the amplitudes of the pion photoproduction process for the 
four physical channels consist as follows 



(mr+|F| 7 p> = -y/2 [if } + F^] , 
(pir~\ Fi 1 771) = [F® - F«] , 

{pAf\ip) = f^+fP, 

(wr°| Fi | 7 n) = F. (+) -F 4 {0) . (B.16) 

The amplitudes F^ 2 ' and F^ 2 ' referring to the final irN isospin states with isospin 
I and ~, respectively, are defined by 

F im = F W_ F H (1 = f)f 

F (i/2) = F ( +) + 2F (-) ( I = iy ( B17 ) 

The isoscalar amplitudes F^ leads only to a 7riV states with isospin |. In terms of 
these, one can define the neutron and the proton multipole amplitudes with total 
isospin I = I as follow 

p(l/2) „(0) _ l F (l/2) 

i,n i 3 ' 

^(1/2) = F (0) + 1 F (V 2)) (B . 18) 

where the subscript p (n) denotes a proton (neutron) target. The four charge chan- 
nels in turn are given by 

F^-n^) = V2[F^-\F^\, 

F,(7n-pO = v/2 [F<f > + iff /2) ] , 

fite-pO = F^ + ^Ff/ 2 ), 

Fiip/n — > n7r°) = - F^ + 2 -F?' 2) , (B.19) 

where the isospin I = | multipole amplitudes, F^ 2 , are the same for both protons 
and neutrons. 

The partial wave decomposition of defines the multipole amplitudes A i± (W), 



B.2 Multipole Decomposition of Amplitudes 



S3 



where i± is the 7riV-orbital angular momentum £ and total angular momentum 
{£ ± |) as well as A = E (A = M) may characterize an electrical (magnetic) 
transition. The '+', '-' determine whether the nucleon spin | must be added to or 
subtracted from the orbital momentum to give the total final state momentum. One 
finds the following relations between the multipole and CGLN-amplitudes ES] 



E, 



+1 



2£ + 2 J_ x 

e+i 



dx 



P e {x)F 1 - P i+ i(x)F 2 + ——(P^x) - P l+x {x))F c 



2£+l 



2£ + 3 



(P e (x) - P t+2 {x))F A 



-i 



uL dx 

4 



2£- 1 



2£ + 2 



P e (x)F 1 -P i - 1 (x)F 2 + 



(P t {x) - P t .. 2 (x))F4 



£+1 
2£ + l 



{P e+ i(x) - P e -t{x))F 3 



+i 



dx 



P t (x)F 1 -P t+1 (x)F 2 + 



— dx 

2£ 



-P i (x)F 1 + P i „ 1 (x)F 2 + 



2£+l 
1 



(P e+1 (x) - P e ^(x))F 3 



2£ + l 



{Pi-x{x) - P e+ i(x))F 3 



;b.2o) 



Here the Pe(x) are Legendre polynomials of the first kind. 



The partial wave decomposition of the CGLN-amplitudes into multipole amplitudes 
Eg±, Me± corresponding to good parity and angular momentum states is very con- 
venient for the description of the resonance excitation process: 

Ft = ^{(£m £+ + e £+ )p; +1 + [(£ + i)m £ _ + ^_]p;_ 1 } ) 

1=0 

oo 

F 2 = ^[(i'+l)M £+ + £M,_]P;, 

1=0 

oo 

1=0 

oo 

F 4 = Y d {M t+ -E t+ -Mt.-E l .)P'^ (B.21) 

1=0 

where the P[ and P" are the derivatives of Legendre polynomials. 
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Appendix C 
Field Operators 



This appendix is concerned with the field operators for nucleon, pion and photon. 
For more information about field operators and the derivation of commutator and 
anti-commutator relations we refer to |51j . Quantized nucleon, pion and photon 
fields are used in the determination of the vertices in chapter El The field operator 
of the nucleon has the form 



with the following anti-commutator relations for the creation and annihilation op- 
erators, respectively 




(C.l) 



or 




(C.2) 



{d(p 1 ,m s ,m t ),d}(p 2 



{b( y p 1 ,m s ,m t ),b^(p 2 ,m s ^m t >)} = (2%) 3 5 3 (p 2 - p[)S, 



? 2 ,m s >,m t >)} = (2tt) 3 (5 3 (p1 -p[)S, 





(C.3) 



(C.4) 



* M N 



The pion field is given by 




f d 3 q 1 



(C.5) 



J (2vr)3 2o; g - 
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with the commutator relation 

K(ql), al,{q 2 )\ = (27r) 3 5 3 (q 2 - q[)8^ 2u ? (C.6) 

— * 

The photon field in Coulomb gauge k ■ e = is given by 

= E / 70 |r + <$e-**) , (C7) 

with the commutator relation 

al^h)] = (2n) 3 5 3 (k 2 - h)5 m ^ 2cu, . (C.8) 

The effect of creation and annihilation operators on the vacuum state is given by 

tf(p,m s ,mt)\0) = \N;pm s m t ) , b(p, m s , m t )\0) = , (C.9) 
at( g -)|0) = |vr;g» , a„(g)|0) = 0, (CIO) 

aJn ,(fc)|0) = | 7 ;fcm 7 > , am7 (fc)|0)=0. (C.ll) 



Appendix D 



Parametrization of the Deuteron 
Wave Functions 



For practical purposes it is more convenient to work with an analytical parametriza- 
tion of the S- and .D-waves in momentum space, Uo(p) and U2{p), respectively, rather 
than with the discrete sets u (pi) and ^(Pi)- The ansatz for the analytic version of 
the momentum space wave functions of the non-relativistic version OBEPR of the 
Bonn OBE-potential jlH] is 



In Table ID. II we list the coefficients Cj and Di of the parametrized deuteron wave 
functions Uo(p) and v,2(p) of the Bonn OBE-potential. In order to make the Uo(p) 
and v,2(p) fulfill the necessary boundary conditions, the last value of Cj and the last 
three values of Di must be calculated through the following relations [15] 




(D.l) 



(D.2) 



with the normalization 




(D.3) 



n u -l 



c, 



1=1 



D 



-2 



{ml 



™l w -2)( m2 n 




-1 
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Table D.l: Coefficients of the parametrized deuteron wave functions for the Bonn OBE- 
potential (full model) jlSJ. The last Q and the last three Di are to be computed from Eq. 



El K 



n,, 



11 



Q (fm 



-l/2\ 



A (fm- 1/2 ) 



0.90457337 +00 
-0.35058661 +00 
-0.17635927 +00 
-0.10418261 +02 

0.45089439 +02 
-0.14861947 +03 

0.31779642 +03 
-0.37496518 +03 

0.22560032 +03 
-0.54858290 +02 
see Eq. (jTOjl 



0.24133026 -01 
-0.64430531 +00 
0.51093352 +00 
-0.54419065 +01 
0.15872034 +02 
-0.14742981 +02 
0.44956539 +01 
-0.71152863 -01 
see Eq. (ITX!) 
see Eq. (ITX!) 
see Eq. (flT4jl 



x 



-m 



E— + 

l=l l 



TLyj 3 



1=1 



n m — 3 



i=l 



(DA) 



and two other relations obtained by circular permutation of (n w — 2), (n w — 1) and 
n w . The masses are chosen to be 



rrii 



a + (i — l)m . 



with a = y/EfyM^ = 0.231609 fm 1 , where Ej, is the deuteron binding energy, and 
m = 0.9 fm -1 . This choice ensures the correct asymptotic behaviour. In analogy 



with Ref. we use n u 



11. 



Appendix E 



Transformation of Differential 
Cross Section 



In Ref. [22] the differential cross sections for the reaction ^(7, n )np are given in 
the so-called photon- nucleon center-of-mass frame 1 . This frame corresponds to an 
assumption that both nucleons in the deuteron have the same momenta — | in the 
7<i c. m. frame. This means that both nucleons in the deuteron are at rest in the 
laboratory frame of the deuteron. 

In order to compare our theoretical predictions for differential cross sections of the 
reaction d(7, ir°)np with the experimental data from Ref. [22] we need therefore to 
transform the differential cross sections from the ~fd rest frame to the c. m. frame of 
the incident photon and a nucleon at rest. Therefore, we give in this appendix the 
Jacoby determinant which we need for this purpose. 

Since the azimuthal angle remains unchanged by the transformation between the 
laboratory and the c. m. systems, we have 



da da dZ 

(E.l) 



dQ cm dfl lab dZ* ' 

where Z*(= cos 9*) and Z(= cos6 l ® b ) are the pion angle in the cm. and the 7c? 
laboratory frame, respectively. The derivative of Z with respect to Z* gives direct 
the requested Jacoby determinant. 



Using the Mandelstam variable t from Eq. ()B.4jl one can express Z through Z* as 
follows 



<7*. * * 771*, ,* 1 rplab, .lab 

Z = ^~ ^ h ^ 7 ■ (E.2) 

, dab n lab v ' 

^7 <?7r 



1 Throughout this appendix parameters in this frame are labelled by (*). 
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which gives the following expression for the Jacoby determinant 




and pion energy (Eq. IB.10)) in the cm. frame, respectively; and u£ ab , q l £ b and E l ® h 
are the photon energy, pion momentum and pion energy in the laboratory frame, 
respectively. 

Using the Mandelstam variable u from Eq. (|B.4|) . one can express the pion energy 
in the laboratory frame as follows 



where E* N = y + lu* 2 is the nucleon energy is the cm. frame. Using the well 

known formulae for momentum and energies from appendix IB. II (see also 01]) the 
pion energy in the laborartory frame is then given by 



El 



dab 



E* N E% + u;q*Z* 
M N 



(E.4) 




El 



7T 



dab 



AB + ufZC 
A 2 - uol ab2 Z 2 ' 



(E.5) 



where 



A 



(E.6) 




(E.7) 



and 




(E.8) 



Appendix F 
Two-Body Subsystems 



In this appendix we give some details of the NN and irN potential models which 
we use to study the NN and irN interaction in the NN and tcN subsystems. In 
the calculations of this work we solve the three dimensional integral equation of 
Lippmann-Schwinger [Uj to calculate the T-matrix elements of the two-body scat- 
tering and then show how phase shifts are obtained. 

The use of separable potentials in the two-body interaction has greatly stimulated 
the theoretical investigation of the three-body problem in the last three decades. 
Most of the important properties of the three-body system are already well described 
by the simplest type of separable models. These separable models are most widely 
used in case of the ttNN system (see for example jHEj and references therein). For 
the NN and ttN interactions in the iViV- and 7riV-subsystems we use a specific 
class of separable potentials [T5J EH1 EH] which historically have played and still play 
a major role in the development of few-body physics and also fit the phase shift 
data for both NN- and 7riV-scattering. The formalism to obtain the full two-body 
T-matrix and the scattering phase shifts is given in this appendix in more detail. 
Momentum space separable potentials for NN- and 7riV-interactions are also given. 



F.l The NN Subsystem 
F.l.l NN Scattering Equation 

In this section we discuss how the T-matrix for iViV-scattering which we use as input 
in our predictions in chapter El is calculated. Two-nucleon scattering is described by 
the Lippmann-Schwinger equation which given graphically in Fig. IF. 11 



In partial wave decomposition, represented by the orbital angular momentum £, 
total angular momentum J, spin s and isospin t, the T-matrix of the two-nucleon 
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^1 = 51 ♦ (v) S (0 

Figure F.l: The graphical diagram of the Lippmann-Schwinger equation for NN inter- 
action. 

scattering may be written as 

{i'J>s't'\T{p',p-E)\Ust) = 8 J ,j5 s r s 8 vt % J e f\p\p-E), (F.l) 

where p and p' are the relative three-momenta of the two interacting nucleons in the 
initial and final state, respectively; and p=|p| and p'=\p'\. E denotes the energy of 

the two interacting nucleons in the c. m. frame and is given by E = j^-, where po is 
the magnitude of the initial relative momentum (c. m. on-shell momentum) which 
is related to the laboratory energy by 

Tiab = ■ (F.2) 

M N 

For a given NN potential V, the T-matrix is obtained from the Lippmann-Schwinger 
equation 02] which in partial wave decomposition reads 1 

%(p',p;E) = V e (p',p)+ dkk 2 V e (p',k) N Mk,p;E) (F.3) 

Jo M N E - k 2 + ie 

for single channels. For the coupled channels it reads 
%e'(p',F,E) = V U '(p',p) 

+ dkk 2 V u „(p', k) N %„ e (k, P] E) , (F.4) 

~ Jo M N E - k 2 + te 

where k is the relative three-momenta of the two interacting nucleons in the inter- 
mediate state and k=\k\. 

Using the identity 
1 



- = V in5(x-Xo), (F.5) 

X — Xq -\- Z6 X — Xq 

where V denotes the principle value integral, one gets the T-matrix elements in 
partial wave decomposition 

f°° M 
%(p', P] E) = V e (p',p)+V dkk 2 Ve(p',k)^-%{k,p;E) 

Jo Po ~ K 

-li7rM N poV i { P ',po)%i P o,p; E) (F.6) 



1 We drop the indices Jst for simplicity. 
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for single channels and 



T w (p',p;E) = V u >(p',p)+ V dkk 2 V u „(p',k)^^%„ e ,(k,p;E) 

i"=j±i Jo Po fc 



■-iTrM N poV U '(p',po)%£'(po,P] E) (F.7) 



for coupled channels. These one-dimensional integral equations, Eqs. (|F.6|) and 
(jF.7|) . can be solved numerically by means of a matrix inversion algorithm. 



Let us first consider the partial wave one-dimensional integral equation for uncoupled 
channels given by Eq. (jF.6|) . The integral should now be approximated by a n-point 
integration routine 

n 

dkF{k) = J^Fik^Wi, (F.8) 

o i=l 

where ki and u>, are the Gaussian integration points and weights, respectively. Since 
the original Gauss points, /Q (with weights Wj) are in the interval (—1,-1-1), we have 
to map them to the interval (0,oo), the range of our integration. For this purpose 
we use the mapping 

h = B tan j^(/Q + l)} , (F.9) 
with the new weights 



4 Icos^K^ + l)] 

where B = 400 MeV is used. For the Gauss integral we use n — 16 mesh points 



To solve Eqs. (|F.6|) and (|F.7j) . the principle value integration has to be replaced by 
a smooth integrand, which can be achieved by adding a term of measure zero 

T t (p',p;E) = V t (p',p)+ dk^-{k 2 V e (p',k)%(k,p;E) 

Jo Po ~ K 

-plVe(p',Po)%(p ,p;E)} - -inM N p V E (p',p )T e (po,p; E). (F.ll) 

All of the n integration points, ki, k<z, , k n , are required to be unequal to po- If 

we call po the (n + 1) point (po = Pn+i), then for p = p Eq. f)F. 1 1 j> can be rewritten 

as 

V e (pi,p n+1 ) = T e (pi,p n+1 ; E)-^2 2 MiV ^ fc 2 {p 2 k Vi(pi,p k )%(p k ,p n+1 ; E) 

k=l Pn+l Pk 

-p 2 n+1 V e (pi,p n+l )%(p n+1 ,p n+1 ;E)} 

+-inM N p n+ iV e (pi,p n+1 )%(p n+ i,p n+1 ; E) . (F.12) 
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Thus, the scattering equation can be written in matrix form as 

n+l 

Ve(Pi,p n +i) = ^Aiip^p^Teipj^n+^E) , (F.13) 

3=1 

where the dimensionless Ae(pi,pj) matrix is given by 

MjyWjp 2 : 
Pn+l Pj 

Wk \ nt 



+ { ( XI ~2 ) + T f ^ArPn+lV^(pi,p n+ l)5 n+ ij 



(F.14) 



Choosing the momentum grid such that pj ^ p n +i V j < n, the matrix A(pi,pj) is 
nonsingular and can be inverted to obtain the on-shell or half-off-shell T-matrix 2 

71+1 

,Pn+l)- ( F - 15 ) 

i=i 

The extension to coupled channels is straightforward. One simply combines the 
points to form a larger (2n + 2) x (2n + 2) dimensional matrix 

2n+2 

V(Pi,Pj) = ^A{p i ,p m )%{p m ,p j ;E). (F.16) 

Here the i, j and m labels include both the points k±, k2, , k n+ \ and the I- value. 

For example, for I = J — 1 we take the 1 < i < (n + 1) points as ki, k.2, , k n+ \. 

For t = J + 1, the label % ranges as (n + 2) < i < (2n + 2) with identical k- values, 

k\i k2, , k n+ \ = k n+ 2i k n+ 3, , /C2n+2- 

So far, we have never mentioned the total isospin of the two-nucleon system, t (which 
is either or 1). The reason for this is simply that t is not an independent quantum 
number. That is, owing to the antisymmetry of the two-fermion state, the quantum 
numbers £, s, and t have to fulfill the condition 

(-iy+ s+t = -i. (F.i7) 

Thus, for a given £ and s, the total isospin t of the two-nucleon state is fixed. 



2 For the matrix T(a, 6; E) we mean by "on-shell" that we are considering elastic scattering. This 
means that the matrix is "on-shell" when E a = Ei> — E, where E is the non-rclativistic kinetic 
energy. If E a ^ E b = E then the matrix is " half-ofT-shell" and it is " full-off-shell" if E a ^ E b ^ E. 



F.l The NN Subsystem 



95 



F.1.2 Separable NN Potential Model 

Now we outline the model which we use to study the NN interaction in the NN- 
subsystem. The EST method jHZj for constructing separable representations of mod- 
ern NN potentials has been applied by the Graz group |IB1 EHj to cast the Paris 
potential [HE] in separable form. In the meantime these so-called PEST potentials 
have become of great use in introducing the features of advanced meson-exchange 
theory jHH] into calculations of few-body systems like, e.g., of nucleon-deuteron scat- 
tering [20] • They have the general form 

N 

V u >(p',p) = J2 9dP') K 9i>M > ( F -!8) 

where iV and A specify the rank and the strength of the separable interaction, 
respectively, in each partial wave. These separable interactions represent a good 
approximation of the on-shell as well as off-shell properties of the original Paris 
potential and show a good fit to the modern NN data base. 

The form factors of these separable potentials for each partial wave are given by 

= E (^+^ n )^-i ' (2<t<3) (F.19) 



4 Cf ^+ 2 (™- 1 ) 

9 " (p) = S^+«j' + "- 2 ' (e - i] 

where C and (3 are fit parameters of the separable interaction. The values of these 
parameters in each partial wave are given appendixO Let us now discuss the partial 
wave states in detail. 

In case of the 1 Sq partial wave the Graz group j3H| have constructed a modified 
PEST potential of rank-3. The parameters of this potential are given in Table lG~T1 
for both the (p-p) and (n-p) systems. Fig. IF.2I shows that the phase shifts of the 
(p-p) system 3 are the same for both PEST3 and the original Paris potential. 

For the 3 Pq, 3 Pi and l P\ partial waves, they constructed a rank-2 separable ap- 
proximations. The potential parameters for these three partial waves are given in 

3 The Paris potential is given only for the (p-p) system in the 1 Sq partial wave. 
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Tiab [MeV] T lab [MeV] 

Figure F.2: 1 Sq phase shift using a rank-3 separable potential from 49 (solid curve) in 
comparison with the original Paris potential (dotted curve) which is given only for the 
(p-p) system. The data points are from SAID [2S] (Solution: SPOO). 



Table In this Table we give the modified parameters of the 3 Pq partial wave 
from Ref. jlHj which we use in this work. The on-shell properties of the PEST 
parametrization for 3 Pq are evident from the phase shift for 3 Po in the left panel of 
Fig. IF .31 Up to Ti a b ~ 70 MeV the approximation to the Paris potential is rather 
good. If we go to higher energies, there is some discrepancy in the phase shifts 
around Ti ao « 100 — 180 MeV, but the situation improves again at higher energies. 
For the 3 Pi and 1 Pi partial waves the separable potential is good up to Ti ab < 200 
MeV. For T lah > 200 MeV the deviation in the 3 Pi and 1 P i partial waves from the 
original Paris potential increases strongly. 

For the 1 D 2 and 3 D 2 partial waves the Graz group constructed a rank-2 PEST ap- 
proximation. The properties of these potentials are practically the same as discussed 
above for the uncoupled P waves. The on-shell behaviour of the potentials, whose 
parameters are given in Table IG.3| are evident from the phase shifts shown in Fig. 

E2D 

The coupled 3 Si- 3 D 1 partial wave state which acting in the deuteron is crucial for a 
proper treatment of the NN interaction. But untill now, no satisfactory description 
of all its aspects has been achieved by means of separable forces. The Graz group 
tried to remedy this situation by constructing a rank-4 PEST approximation to the 
Paris potential. This separable interaction, whose parameters are given in Table RT4"! 
yields a correct reproduction of the deuteron properties as well as a fair description 
of the scattering domain as shown in Fig IF. 51 
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100 200 300 100 200 300 100 200 300 

Tiab [MeV] T lab [MeV] T lab [MeV] 

Figure F.3: Same notation as in Fig. lF.2| but for the s Pq, l P\ and 3 Pi and single channels 
using a rank-2 separable potential. 



For the coupled 3 p2- 3 p2 partial wave, they constructed a rank-3 PEST parametriza- 
tion, whose parameters are given in Table IG.51 The phase shifts and the mixing 
parameter are shown in Fig. IF.6I On the other hand, they designed a rank-4 PEST 
parametrization for the coupled 3 D 3 - 3 G3 partial wave. The parameters of the PEST 
potential are given in Table RT6l The on-shell properties are shown in Fig. IF. 71 they 
exhibit a good agreement with the Paris potential. 

From the foregoing presentation of the PEST interactions it is clear that these 
separable representation of the Paris NN potential are likely to be relevant in few- 
body problems. 



F.1.3 NN Phase Shifts 

The full on-shell S"-matrix in partial wave decomposition is related to the on-shell 
partial wave T-matrix by 

S(E) = l + 2iT(E), (F.20) 

with 

T(E) := ~7r Po M N T{E), (F.21) 
in the nonrelativistic kinematics which we use in this work. 

The scattering phase shifts for a given partial wave can be calculated from the on- 
energy-shell i^-matrix (a 2x2 matrix) (ZUIZ21 which is expressed in terms of the S- 
and T-matrix by 

i _ Q T 

K := i- = — . (F.22) 

1+S 1+iT K ' 
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Figure F.4: Same notation as in Fig. IF.2l but for the 1 L>2 and 3 -D2 single channels using 
a rank-2 separable potential. 



For the single channels, the scattering phase shifts are then determined by 

tan 5 = $leK. (F.23) 

For the coupled channels, a unitary transformation is needed to diagonalize the 
two-by-two coupled ^-matrix. This requires an additional parameter, known as the 
mixing parameter e. This mixing parameter and the phase shifts of the coupled 
channels are given in terms of the .fT-matrix by 



sin(2e) 



(2^eK ) 2 



[1 + ($leK ) 2 - MT+KefC] 2 + i^eK + + $leK_) 2 



(F.24) 



and 



25± 


= arcsin ^ 


$leK + + $leK_ 
2^eK 


where 






K 


:= K J+1 ,j. 


-l = Kj-i t j + i 


K + 


'■= Kj+i,j- 


hi 5 






-1 ■ 



sin(2e) ± arcsin 



( 2^eK 



tan(2e" 



(F.25) 



(F.26) 



In this work, all phase shifts are in the Stapp-Ypsilantis-Metropolis (SYM) or bar 
convention |73j . 
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100 200 300 100 200 300 100 200 300 

Tiab [MeV] T lab [MeV] T lab [MeV] 

Figure F.5: Same notation as in Fig. IF. 21 but for the coupled channel 3 Si- 3 -Di using a 
rank-4 separable potential. 

F.2 The nN Subsystem 

Analogous to the case of iViV-scattering, the formalism to obtain the full T-matrix 
and the scattering phase shifts for 7riV-scattering is outlined in this section, but not 
in more details since we give enough details in case of the iViV-scattering. A partial 
wave decomposition of the separable irN potential is performed in this section. The 
nN scattering equation is also given. 



F.2.1 7riV-Scattering Equation 

For a given tcN potential V, the T-matrix for nN scattering is obtained from the 
Lippmann-Schwinger equation which in partial wave decomposition, specified by the 
orbital angular momentum £, total angular momentum J and isospin t, reads 4 

POO 

T(p',p;E) = V(p',p)+ dkk 2 V(p',k)g nN (kiE)T(k,p;E), (F.27) 

J o 

with the ttN propagator 

Mk,E) = = _ 7TT ~ J-, 7TT . , (F.28) 

E — E^^k) — E N (k) + le 

where p, k and p' are the irN relative momentum in the initial, intermediate and 
final state, respectively; andp=|p|, k=\k\ and p'=\p'\. E = a/ + p\ + a/ m\ + p^ 



4 The indicies £Jt are also dropped here for simplicity. 
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Figure F.6: Same notation as in Fig. IF .21 but for the coupled channel 3 P2- 3 -^2 using a 
rank-3 separable potential. 



denotes the total collision energy with on-shell momentum p in the c. m. frame. 
For our calculations we use relativistic kinematics for both pion and nucleon, thus 

K(k) = V™1 + k 2 , 



E N (k) = ^M 2 N + k 2 . (F.29) 
The partial wave T-matrix for ttN scattering is then given by 

T{p',p;E) = V(p',p)+ ^ 'dkk 2 V(p\k) 2 T(k,p;E), (F.30) 

Jo p — k + is 

where 

r(k] [E N ( Po )+E N (k)} [E n ( Po ) + E n (k)} 

1 1 [E N (po) + E N (k) + E„(p ) + E v (k)] • 1 ' ] 

To transform this equation into a principle value integral equation, one uses the 
identity 

' V = 1 2 - mb(x 2 - xl) , (F.32) 



•Xj n T" v C Jo Ju . 



to get the matrix elements in partial wave decomposition as 

T(p>,p;E) = V(p',p)+V [°° dkk 2 V(p',k)-^^T(k,p;E) 

Jo Po ~ fc 



-^•7rp G(po)^(p , ,Po)r(po,P; J B) • (F.33) 
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Figure F.7: Same notation as in Fig. IF, 21 but for the coupled channel 3 Z?3- 3 G3 using a 
rank-4 separable potential. 



This one-dimensional integral equation can be solved numerically for a given ttN 
potential model V(p',p) by using the matrix inversion method which is explaind in 
details in case of iViV-scattering in the previous section. In this work we use the 
separable nN potential model of Nozawa et al. [15J. This model is given in more 
details in the forthcoming section. 



F.2. 2 Separable ttN Potential Model 

In the case of the ttN scattering a large number of dynamical models have been 
developed over the past few years (see for example [13 EH ESI)- Most begin with a 
separable potential which is iterated in a Lippmann-Schwinger equation to give the 
scattering amplitude, from which phase shifts and observables are obtained. In this 
work we use the model of Nozawa et al. ^Hj in order to study the ttN interaction in 
the 7riV-subsystem. This model is consistent with the existing unitary description 
of the irNN system and treats the ttN interaction dynamically, with all S-, P- and 
.D-wave ttN phase shifts being well reproduced below 500 MeV. 

For partial wave specified by quantum numbers £Jt, the ttN scattering equation 
takes the form given in Eq. (JE2ZJ). The vriV potential V(p',p) is given in Ref. |15j 
due to both a vertex interaction f° N B and a "background" potential V (which is 
assumed to be of a separable form) by 

V(p',p;E) = J2 £n,b9b(E) f°B,„ N + V, (F.34) 

B=N,A 
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where g%{E) is the free propagator in subspace 5 B and takes the form 

9°b = ~p "~ 5 (F.35) 
with m B is the bare mass of the baryon B. 

For P\i and P33 channels, both terms in the right-hand-side of Eq. (|F.34J) are taken 
into account. Thus the total potential matrix element can be written in terms of 
baryon pole and non-pole parts as follows 

V(p',p;E) = fo(p')g B (E)fo(p) + ho(p')\oh (p), (F.36) 

where Ao is a phase parameter and it is given in Table IHTTI The form factors ho(k) 
and fo(k) are parametrized as follows 

n umi 



a 2 k m2 
(k 2 + bl) n * 



Mk) = , ^ k £ - (F.38) 



The values of the parameters ni, n,2, ra%, m>2, ai, a<z, b\ and 62 are given in Table 
IH.ll for each partial wave. 



By inserting Eq. (|F.36j) into Eq. (|F.27|) . the irN amplitude can then be written as 

T(p',p;E) = T NP (p',p;E) + f(p')g B (E)f(p). (F.39) 

The non-pole T^-matrix is given by 

T NP (p',p;E) = h (p') T (E) h (p) , (F.40) 

with 

To{E) = i-XoJdkk* \h (k)\ig wN (k,E)- (R41) 

The pole term in Eq. (|F.39J) consists of a dressed form factor f(k a ) (a = i, f) defined 
by 

f(Pa) = fo(Pa) + r (E) h (pa) [ dkk 2 h (k) f (k) Q 7rN (k,E) , (F.42) 



'Subspace 5 = J2b=n a n ^ ® & describes irN scattering without coupling to photons. 
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and a dressed propagator Qb(E) defined by 

9b{E) = E-moB-f dkk* f(k) f Q (k) GMk, E) • (R43) 

For channels other than Pu and P33, the vertex interaction of Eq. (|F.34j) does not 
contribute. The ttN potential is then assumed to be of rank 2 separable form 

V(p',p) = h^p') Ai h(p) + h 2 (p') A 2 h 2 (p) . (F.44) 

The form factors are parametrized as 



a 2 k m2 



As before, the parameters n%, n 2 , mi, m 2 , a±, a 2 , h\ and b 2 are given in Table ITTT1 
for each partial wave. 

By inserting Eq. (|F.44J) into Eq. (|F.34jl . one obtains the following analytic solution 
for the T-matrix 

T(p',nE) = h 1 (p')m(E)hi(p)+hi(p / )r 12 (E)h 2 (p) 

+h 2 (p') r 21 (E) h^p) + h 2 (p') r 22 (E) h 2 (p) , (F.47) 

where 

me) = n : TT :;r :z7 ma- . ( R48 ) 



X 1 (l-X 2 H 2 ) 


(1-X 1 H 1 )(1-X 2 H 2 )- 


X\X 2 E 2 2 


t 2 i(E) 




XiX 2 H\ 2 




(1 - X^l - X 2 H 2 ) - 


X\X 2 Hl 2 


Aa(l-Aifl-i) 




(1 - X 1 H 1 )(1 - X 2 H 2 ) - 


X\X 2 Hf 2 



(F.49) 

me) = - — x rr ;r7'^ i; , , ^2 > ( F - 5 °) 

with 

E x = I dkk' 2 Q %N (k,E), (F.51) 



tf 2 = J dkk 2 \h 2 (k)\ 2 g %N (k, E) , (F.52) 
#12 = J dkk 2 h{k) h 2 {k) GMk,E). (F.53) 
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Tiab [MeV] T lab [MeV] 

Figure F.8: The ttN phase shifts of the S partial waves obtained from the LS equation 
using the separable potential model of Nozawa et al. |15j shown versus the pion laboratory 
energy T/ Q & in MeV. The data points are from the VPI partial wave analysis [2B] (Solution: 
SM99). 



F.2.3 ttN Phase Shifts 

Once the scattering equtaion given in Eq. (|F.33J) is solved for the on-shell T-matrix, 
one can obtain the scattering phase shifts. The information about the scattering 
process is commonly represented by phase shifts 5(E) with the on-shell S'-matrix in 
each channel denned by 

S(E) = e 2t5{E \ (F.54) 



or 



S(E) = 1 + 2i e i5{E) sin(5(E)) , (F.55) 

where E(p) = a/ Mjy + p 2 + a/ rr^ + p 2 is the invariant total energy of the two 
interacting praticles. The partial wave on-shell S'-matrix is related to the partial 
wave on-shell T-matrix by 



S(E) = 1 - 1m p(p) T(E) , (F.56) 
where the density of states p(p) is given by 
( \ J? pE N (p)E n (p) 

dp 
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Figure F.9: Same as Fig. IF.8l but for the P partial waves. 



Combining Eqs. (jF.54|) and (|F.56|) gives 

-2tt pip) dte(T(E)) 



tan(25(E)) 



(F.58) 



1 + 2vr p{p) %m(T(E)) ' 
Alternatively, we can express the on-shell T-matrix in terms of the phase shifts as 



T(E) 



whence 



-1 



e i5{E) nn(6(E)), 



7T 



Pip) 



tan(5(£)) 



3m{T(E)) 
®.e{T{E)) ' 



(F.59) 



(F.60) 
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Figure F.10: Same as Fig. IF. 81 but for the D partial waves. 



The phase shifts calculated from the dynamical model of Nozawa et al. [T^j for the 
more important ttN partial waves are shown in Figs. lF.8llF.10l with the corresponding 
parameters being given in Table |H~T1 We see that the Lippmann-Schwinger equation 
gives a good description of nN phase shifts below 500 MeV. We note in particular 
the perfect resonance shape of the P33 phase shift (see Fig. IF.Qjl corresponding to 
the A(1232) resonance of the nucleon (5=90° at 7} i, ~ 180 MeV). The steep rise 
of the Pu phase shift for — > 500 MeV is an indication of the Roper resonance 
N(1440). In addition there is some background scattering due to interactions in the 
Su and 5*31 partial waves. 



Appendix G 

Parameters of the TV TV Separable 
Potential 

We give in this appendix the parameters which we use for the separable NN inter- 
action of the Graz group jlHl EH] in each partial wave. 
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Table G.l: Parameters of the PEST3 potential in the 1 Sq partial wave. 



(3 (fm" 1 ) C (fin ) A (MeV fm- 1 ) 



Piv 


= 1.1115753 


Cn= 


-572.29230 


An- 


-0.001021062 




=2.0212320 


Cl2=- 


-313.62360 


Al2 = 


0.022965747 




=2.6434278 


Cl3 = 


12814.143 


Al3=' 


-0.013907568 


Pu- 


=4.0877911 


Cl4— ■ 


-28477.672 






P21- 


= 1.0297944 


C 2 i= 


-18.494719 


A 2 i= 


0.022965747 


&2 = 


= 1.5435994 


^22=- 


-53.628348 


A 2 2=' 


-0.805816820 


(323- 


=2.6129194 


C*23 = 


925.42408 


A 2 3 = 


0.825835300 




=4.0837929 


C 24 =- 


-2040.0020 






Psv 


=0.8995186 


C*31 = - 


-2.8837445 


A31— 


-0.013907568 


A»2 = 


=2.6843334 


C32=- 


-84.344154 


A32 = 


0.825835300 


(333 


=3.1505869 


C33— 


1163.1547 


A33=' 


-1.295713900 


Pu- 


=3.9826288 


C34 — 


-2018.0523 






Piv 


= 1.1115753 


Cn= 


-773.80000 


An- 


-0.000528400 


Pl2 


=2.0212320 


C12— ■ 


-424.05244 


Al2 = 


0.015471054 


013 = 


=2.6434278 


C*13 = 


17326.082 


Al3=' 


-0.008800256 


Pw 


=4.0877911 


Cl4=- 


-38504.839 






P2V 


= 1.0297944 


C2I— - 


-18.494719 


A 2 i= 


0.015471054 


P22- 


= 1.5435994 


C*22=- 


-53.628348 


A22— 


-0.742248760 


P23- 


=2.6129194 


C23 = 


925.42408 


A 2 3 = 


0.770689240 


#24= 


=4.0837929 


C24=- 


-2040.0020 






Psv 


=0.8995186 


C3I— - 


-2.8837445 


A31— 


-0.008800256 


P32- 


=2.6843334 


C*32=- 


-84.344154 


A32 = 


0.770689240 


P33 


=3.1505869 


^33 = 


1163.1547 


A33=' 


-1.252847400 


Pu~ 


=3.9826288 


C34 — 


-2018.0523 







109 



Table G.2: Parameters of the PEST2 potential in the uncoupled P waves. 



(3 (fm" 1 ) C (fin ) A (MeV fin" 3 ) 





l. 8691379 


C u =- 


-189.74407 


An= 


r\ -1 y( AA f A AA 

-0.14665888 


Pl2 — 


3 n4QQ8fifi 




9nQfi 7981 


Al2 = 


fl 40.^1 ^78^ 


M13 — 


3 ^8^7^ 

O.tJtJOtJ 1 tJZ, 


n,„— 


^flfifi 9fiQ4 






Pl4 — 


J. . ±Ocf\J£c/£ 


O14 — 


1 D9 7^84 






/?21 = 


2.5512467 


C21 — 


59.861550 


^21 = 


0.40513783 


#22 = 


2.6651826 


C22=- 


-681.32014 


^22=- 


-5.47893200 




2 9537139 


n — 

l - y 23 


5265 5015 








3 3827130 




( - y 24 — ' 


-6825 7570 






#11 = 


l. 6225099 




76.181266 


An= 


0.66973064 


Pl2 — 


9 471 9881 


( -^12 — 


^81 78^QQ 
-001 . 1 ooyy 


Al2=- 


n ^990.^1 


Pl3— 


9 71 9DD c iQ 


^13— 


1471 78rn 






Pl4 — 


3 1 <S498rn 


014— 


1 fi4fi 41 fi8 






#21= 


1.7697160 


C^21 — 


61.335643 


A2i=- 


-0.52203105 


#22= 


2.4049729 


C22 = 


1348.5020 


A22 = 


0.64075503 


#23= 


1.8631357 


C*23 = - 


-842.10204 






#24 = 


4.1084782 


C24=" 


-1604.2632 








1.5652079 


Cn= 


89.398777 


An= 


0.44514549 


#12= 


2.0058818 


C12— • 


-208.13553 


Al2=- 


-0.34067880 




2.4100051 


Cl3= 


425.34413 






#14= 


4.2981726 


Cl4=- 


-630.35239 






#21= 


1.5463024 


C21 — 


43.477263 


A2i=- 


-0.34067880 


#22= 


1.6196926 


C*22 = 


166.46466 


A22 = 


0.44501514 


#23 = 


2.0004143 


C23=- 


-404.03968 






#24 = 


2.1796197 


C*24 = 


528.13575 
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Table G.3: Parameters of the PEST2 potential in the uncoupled D waves. 



(fm" 1 ) C (fm°) A (MeV fin" 1 ) 





1.1690501 


C n =- 


-7.6541157 


A n = 


-1.90862530 




1.3354235 


Cl2 = 


3.8397789 


Al2 = 


1.00851240 


013= 


3.6311131 


Cl3— ' 


-348.08899 






014= 


2.9623562 


Cl4 = 


330.93549 






021 = 


0.8715792 


C21 = 


-1.0867868 


^21 = 


1.00851240 


022 = 


1.2853211 


C22=' 


-22.731322 


^22=' 


-0.97543827 


023 = 


3.2457085 


C23 = 


265.30215 






024 = 


3.8839614 


C 24 =- 


-303.60967 






011 = 


1.2055706 


C n =- 


-40.322168 


A n = 


-0.26674102 


012 = 


1.2960008 


Cl2=- 


-21.021097 


Al2 = 


0.14305117 


013 = 


1.6239653 


Cl3— 


533.48026 






014 = 


1.4452618 


Cl4=- 


-436.55162 






021 = 


1.5953203 


C21 = 


-245.06021 


A21 = 


0.14305117 


022 = 


0.3626991 


C*22 = 


3.0950339 


A22=' 


-0.15817248 


023 = 


1.7370869 


C*23 = 


299.77237 






024 = 


0.5873583 


C7 2 4= 


42.826073 
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Table G.4: Parameters of the PEST4 potential in the coupled 3 S\- 3 Di partial wave. 



# (fm" 1 ) C (fm°) A (MeV fm- 1 ) 



#n= 


1.6855291 


C u =- 


-51.009539 


An= 


-0.254222420 


#12 = 


3. 9205339 


<^12 — 


765.24390 


A12— 


-0.171787720 


HIS — 


F5 7636840 


C*13=' 


-3363 4996 


Aiq— 


071 72381 4 


filA — 


6 941 969^ 


Cl4— 


2783 6926 


Ai4 = 


00^1 0^)01 1 


#21 = 


1.7877900 


C21=' 


-68.041320 


A2i=- 


-0.171787720 


#22 = 


2.1759210 


C*22 = 


488.86119 


A22 = 


0.087724390 


#23 = 


2.4705717 


C*23=- 


-1199.1640 


A23 = 


0.094518781 




2 7303931 




871 79604 


A24 ' 


-0 079244845 

V 7 • V 7 1 '.7 J- J- v7 J- .J 


#31 = 


1.5971102 


C3I— " 


-89.517302 


A31 = 


0.071723814 


#32 = 


9.9678931 


C*32 = 


2582.2164 


A32= 


0.094518781 


P33 — 


a t^o/i am 1 


7^ 

^33— 


-oU4z.oozo 


A 33 — 


-u.uouo4oyyz 




2.1206347 


^34 


1235 3473 

-L -.7 7 • -.7 J- 1 -.7 


A34 * 


-0 034997609 

V 7 • V 7 " 7 -L '.7 '.7 1 V 7 V 7 '.7 


#41 = 


4.0268195 


^41 = 


-382.75856 


Mi— 


0.005105011 


#42 = 


5.0466356 


642 = 


1005.8091 


X42— 


-0.079244845 


#43 — 


Z.o i yoyol 


L/43 — 


zzyo. f 0o4 


A43— 


-u.Uo4yy /ouy 


0AA = 


2 3953665 


C44 - 


-2718 2533 


A44 


169553430 

V 7 • -L V7 '7 7 .7 '.J i'.J V 7 


#n= 


2.6228398 




-410.15751 


An= 


-0.254222420 


#12 = 


1.8815276 


( - / 12 — 


156.64127 


Al2=' 


-0.171787720 


yl3— 




Cl3= 


888 f)3Q8 c i 
OOO.Uuc/OJ 


^13— 


071 79381 4 


fi,A — 

M14 — 




Cl4— ■ 


-792 1 1 994 


Al4= 


00^1 0^)01 1 


#21 = 


1.9098545 


C21 = 


87.065827 


A21=' 


-0.171787720 


#22 = 


1.1170776 


C*22 = 


20.554605 


A22= 


0.087724390 


#23 = 


1.4259853 


C23=' 


-65.702538 


A23= 


0.094518781 


#24 = 


3.0734784 


C 2 4 — 


-91.839262 


A 2 4=- 


-0.079244845 


#31 = 


2.6431889 


Csi— 


143.62150 


A31 = 


0.071723814 


#32 = 


2.7483825 


C*32 = 


6.4744841 


A32= 


0.094518781 


#33 = 


1.9283993 


C33 = 


-1384.0981 


-^33=- 


-0.030645992 


#34 = 


2.2891433 


^34 = 


1550.1289 


A34— 


-0.034997609 


#41 = 


4.2598369 


Ca\ = 


366.50118 


A41 = 


0.005105011 


#42 = 


2.1463834 


Ca2 = 


263.87202 


A42 — 


-0.079244845 


#43 = 


2.4905282 


C43 — 


-611.66811 


A43 — 


-0.034997609 


#44 = 


2.2930841 


C44 — 


-12.606207 


A44= 


0.169553430 
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Table G.5: Parameters of the PEST3 potential in the coupled 3 P2- 3 -^2 partial wave. 



(3 (fin" 1 ) C (fin ) A (MeV fm" 3 ) 





1.4452936 


di= 


-28.609391 


An= 


0.04672503 


/?12 = 


2.0173835 


Cn — 


-186.15068 


A 12= 


0.41710484 


/5l3 = 


5.4463467 


Cl3= 


-1094.0574 


Ai 3 = 


-0.06092248 


014= 


2.8490316 


Cl4 = 


1022.0712 






#21= 


1.8993563 


c 21 = 


65.794074 


A 2 i= 


0.41710484 


#22 = 


2.4547294 


C*22 = 


173.01447 


A 22 = 


0.25492496 


&3= 


1.5499793 


C2 3 = 


-190.40569 


A 23 =- 


-0.41481703 


#24 = 


5.9773387 


C24 = 


-217.74704 






Au= 


2.9659955 


C 3 i— 


-327.40790 


A 3 i= 


-0.06092248 


# 32 = 


4.9324835 


C 3 2 = 


-2950.6615 


A 3 2=' 


-0.41481703 


# 33 = 


2.9051889 


c 33 = 


4127.4473 


A 33 = 


-0.13593596 


#34= 


2.0463803 


C 3 4 = 


-1298.0249 








1.1667491 


Cn= 


0.96744407 


An= 


0.04672503 


#12= 


1.5406619 


Cl2 — 


166.30912 


Al2 = 


0.41710484 


As= 


2.3802869 


Cl3— 


-119.42176 


Ai 3 =- 


-0.06092248 


#14= 


4.8335309 


Cl4 = 


366.96981 






/?21 = 


1.0955323 


C*21 = 


-0.6386842 


A 2 i= 


0.41710484 


#22 = 


1.6020596 


C*22 = 


-37.538146 


A 2 2 = 


0.25492496 


#23 = 


4.4585234 


C2 3 = 


-316.97599 


A 23 =- 


-0.41481703 


# 24 = 


1.8254989 


C*24 = 


179.59323 






#31= 


2.1084981 


C31 — 


-47.364933 


A 3 i=- 


-0.06092248 


# 32 = 


3.1977404 


c 32 = 


2035.5781 


A 3 2=- 


-0.41481703 


# 33 = 


3.9573349 


c 33 = 


-5925.4169 


A 33 =- 


-0.13593596 


# 34 = 


4.5120514 


C 3 4 = 


5407.1036 
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Table G.6: Parameters of the PEST4 potential in the coupled 3 D^- 3 Gs partial wave. 



# (fm" 1 ) C (fin ) A (MeV fm- 1 ) 



#n= 


1.7275071 


Cn= 


-8.6379311 


An= 


-0.09899224 


#12 = 


2.549331b 


^12 — 


178.84657 


A12— 


-0.39484562 


H13 — 


1 Q1 1 1 9^4 

X . t_/ X X X JUt: 


C*13=' 


-114 21507 


Aiq — 

A 13 — 


0401 3^28 




4 06691 72 


Cl4— ■ 


-86 040765 


Ai4 = 


1 1 343039 


#21 = 


1.9189644 


C21=- 


-116.07858 


A21 = 


-0.39484562 


#22 = 


2.7381690 


C*22 = 


354.48169 


A22 = 


0.28242764 


#23 = 


1.5396020 


C*23 = 


196.12620 


A23 = 


0.40048398 


A- / 24 


1 9159334 

-L • ' ./ _L 7 '.7 -7 "7 1 


C24 ' 


-419 69957 

-L -L '../ • V 7 7 '.7 7 1 


A24 " 


-0 34677241 


#31 = 


3.2277180 


C3I— " 


-110.31148 


A31 = 


0.04013528 


#32 = 


4.8060158 


C*32 = 


483.29454 


A32 = 


0.40048398 


P33 — 


o.iuuouyu 




44Z.oUUoy 


A 33— ' 


n 1 AOOA 1 01 




3 7214195 


^34 ' 


-903 63957 


A34 " 


-0 38000970 

V 7 • -.7 v..- V 7 V 7 V 7 '7 1 V 7 


#41 = 


2.8426813 


Ca\=- 


-364.88710 


A41 = 


0.11343039 


#42 = 


3.3549103 


C42— 


1848.1193 


A42 = 


-0.34677241 


#43 — 


q 0101 om 


L/43 — 


Z00U.004 1 


A43 — 


-U.ooUUUy 


3aa = 


2 7766588 

■ 1 1 V7 V7 \j v7 v7 


C44 ■ 


-3985 0257 


A44 


50560372 


#n= 


2.4000943 


Cn= 


10.076418 


A n = 


-0.09899224 


#12 = 


1 noon /i o ^7 

1.2830437 


✓ ■ 

W2 — 


-90.394797 


^12 — 


-0.39484562 


/J13— 


1 fifi97199 


Cl3= 


983 81 7D7 


^13— 


fl (14(11 3^98 


M14 — 


1 ^81999 


■ — ■ 


-204 ^2799 


Al4 = 


1 1 343039 


#21 = 


3.9194456 


C21 = 


10.924992 


A21=' 


-0.39484562 


#22 = 


1.5719585 


C*22 = 


82.607174 


A22 = 


0.28242764 


#23 = 


1.6553410 


C23=' 


-132.30266 


A23 = 


0.40048398 


#24 = 


2.0963683 


C 24 = 


47.770331 


A 2 4 = 


-0.34677241 


#31 = 


4.8032746 


C3I— - 


-204.05105 


A31 = 


0.04013528 


#32 = 


3.6976328 


C*32=- 


-751.61370 


A32 = 


0.40048398 


#33 = 


2.7214754 


^33 = 


1886.8354 


A33=' 


-0.14224121 


#34 = 


1.9739909 


^34=- 


-954.91540 


A34 — 


-0.38000970 


#41 = 


1.8399157 


C41— - 


-0.2429913 


A 4 i= 


0.11343039 


#42 = 


2.0748484 


C*42 = 


197.58830 


A42 = 


-0.34677241 


#43 = 


2.2831760 


C43— - 


-411.75699 


A43 — 


-0.38000970 


#44 = 


2.6095142 


C 4 4 = 


229.35025 


A44 = 


0.50560372 
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Appendix H 

Parameters of the ttN Separable 
Potential 



Table H.l: Parameters of the ttN separabel potential of Nozawa et al. ^H] for the ttN 
partial waves. These parameters are determined by fitting the phase shift data |7S] up to 
500 MeV pion laboratory kinetic energy. 



L2t,2J 


t 


mi n\ 


a[ 


b\ 


m 2 


"■2 


«2 


°2 


Ao 


Sll 








3 


100.00 


2.598 


2 


2 


4.9520 


2.877 


-1 










2 


3.0850 


1.806 


2 


2 


1.9250 


1.275 


+1 


Pu 


1 


2 


3 


31.623 


2.665 





2 


0.5793 


1.185 


-1 


Pl3 


1 





2 


0.4269 


1.181 


2 


3 


3.9700 


1.721 


+ 1 


Psi 


1 





2 


1.4730 


1.542 


2 


3 


8.0530 


1.861 


+1 


-P33 


1 





2 


2.7700 


1.415 





2 


1.7780 


1.218 


+ 1 


D 13 


2 





2 


1.6390 


2.165 1 


2 


3 


9.3120 


3.263 


-1 


D 15 


2 





2 


0.2172 


1.175 


2 


3 


1.0110 


1.461 


-1 


D 33 


2 





2 


0.1306 


1.128 


2 


3 


1.0810 


1.972 


-1 


D 35 


2 





2 


0.2270 


1.168 


2 


3 


1.1510 


1.780 


+1 


ai is £ 


dven 


in 


units of (fm) 


2ni+mi+£-\ 


-1 











* ci2 is given in units of (fm) 2n 2 +m 2 +^+ 2 f or anc j p 33 p ar tial waves, otherwise 
it is given in units of (fm) _2n2+m2+ ^ +1 . 

§ b\ and 62 are given in units of (fm) -1 . 

' This value is a missprint in Table 2 of Ref. ^Hj- We would like to thank 
Professor S. Nozawa for giving us the correct one. 
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Chapter 1 

Polarization Observables 

o 

(N 

Polarization observables have the promise of opening a new filed in the electromag- 
netic production of pions from protons and nuclear targets. From this point of view 
■ we give in this chapter the formal expressions of all possible polarization observables 

in the 7c? — > ttNN reaction with polarized photons and oriented deuterons in terms 
of the .M-matrix elements. The importance of this process derives from the fact 
', that the deuteron, being the simplest nuclear system, plays a similar fundamental 

role in nuclear physics as the hydrogen atom plays in atomic physics. 

m 

The initial state for photoproduction of 7r-mesons on the deuteron is characterized 
by the incoming photon polarization m 7 and the deuteron spin projection with 
respect to a chosen quantization axis, in our case for convenience the incoming 
photon momentum k. Thus we associate a frame of reference with z-axis in the 
direction of k. 

The final state is described by the two-nucleon relative momentum p r having the 
spherical coordinates 6 r and r with respect to the chosen frame of reference, by the 
total spin s and its spin projection m s on p r and by the pion momentum q having 
the spherical coordinates 9 n and <j> n . The x-z-plane is determined by the outgoing 
pion, i.e. = 0. This means that the .M-matrix is spherical harmonics depending 
on the pion momentum and its direction. 

Choosing the z-axis in the direction of the incoming photon and isolating the az- 
imuthal dependence of the direction of pion momentum, we obtain the following 
general form for the reaction matrix 

Using parity conservation one can show that the .M-matrix elements obey the fol- 
lowing symmetry relation 

/n(*A») _ ( _\s+ m s+m 1 +rn d (r\(tfi) /i n\ 

^s,— m s ,— m 7) — m d \ ) ^ s,m a ,m 1 ,m d ' K^-^J 
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This symmetry relation reduces the number of complex amplitudes from 24 to 12 
independent ones. For their determination one needs 23 real observables since a 
overall phase remains arbitrary. 

For a given .M-matrix one can compute the cross section for arbitrary polarized 
photons and initial deuterons by applying the density matrix formalism similar to 
that given by Arenhovel j2] for deuteron photodisintegration. The most general 
expression for all observables is given by 

= Ps £ (A^:? m , m , )* M^\m imd Pl Wl Pt d m> dOfr , (L3) 

sm s i ,m~y m d 

s 1 m' B t' ' ,m'^m! d 

where p^ m / and p d mdm i denote the density matrices of initial photon polarization 

and deuteron orientation, respectively, and Q is an operator associated with the 
observable which acts in the two-nucleon spin space. The kinematical phase space 
factor p s is given in pQ. For more details about the density matrices we refer to [2]. 

For a given .M-matrix all polarization observables can be expressed in terms of the 
following quantities 



.m E 



v IM = > > (-y- m *i 



l-m', - I 1 1 I 



m' d m d sm s !,m 1 

x / p. (M£„)* -^U^ *V > (L4) 



and 



wim = ± E E m 1 -* ' 



i i / 

m' d m d sm s t,m~f x 

x / P. (MS^)* ^;i, mT _ m - dfift , (1.5) 



where I = y/2I + 1. 

The unpolarized differential cross section is then given by 
da 

V w ■ (1-6) 



d^dq 

The photon asymmetry for linearly polarized photons is given by 
The target asymmetries are given by 



3 



The vector target asymmetry: 



da 



• The tensor target asymmetries: 

da 



T 2M m = [^-j (2 - 8 M0 ) MeV 2M , (M = 0,l,2). (1.9) 

The photon and target asymmetries are given by 
• Circular: 

( d \ _1 

T^ M (e n ) = l^f^J (2 - Smo) $teV 1M , (M = 0,1) 

da 
dQ n dq 

Longitudinal: 

T ' m( ^ } = (d^dq) lWlM ' ( M = ' ±1 )' 

TU0«) = W ™, (M = 0,±l,±2). (1.11) 



T 2m(^) = (-^H 2Smy 2M) (M = 0,l,2), (1.10) 
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